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Various problems in physics, for instance, wave propagation and heat conduction in 
materials with memory effects are of great interests. Various models in continuum 
mechanics describe the mechanical behavior of solids. In elasticity, an instantaneous 
response of material is treated as an adequate approximation. Finite elasticity is 
a part of elasticity theory and does not put any a priori limitations on displace- 
ment and strains. Such physical problems can be modeled as initial boundary value 
problems for partial integrodifierential equations. Our main concern is viscoelastic 
models, where viscoelastic theory may be considered as a theory taking into account 
of the effect of the entire deformation history on stresses at current instant. In gen- 
eral materials exhibit two different kinds of behaviors, liquid like and solid like. Here 
we consider the conservation of momentum and energy in viscoelastic media with 
finite deformations. Such initial boundary value problems may be reformulated as 
abstract integrodifferential equations in abstract spaces such as Hilbert spaces or 
Banach spaces. The main advantage of studying these problems in abstract spaces 
is that we may only have to concentrate on invariant properties of the problems and 
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need not worry about unnecessary details of a specific problem, and the established 
results can be applied for the whole class of problems. 

We consider the abstract formulations of various physical problems of viscoelas- 
tic solid models with short and long memory and related problems. For a material 
with memory, the stress at an instant of time depends in some fashion on the history 
of the strain up to the instant of time. The governing equation of motion of a mate- 
rial with memory can be represented as an integrodifferential equations formulated 
in an abstract space, such as a Hilbert space or a Banach space. By considering 
their abstract formulations we may concentrate more on the main properties which 
remain invariant for the whole class. The rich theory of functional analysis may 
be used to establish wellposedness and the convergence of various type of solutions 
to the problems under consideration. Our main tools from the functional analysis 
are the semigroup of operators in a Banach space, theory of monotone (accretive) 
operators in a Banach space and the method of semidiscritization in time. The main 
importance of the theory of semigroup appears when we handle the problem with 
unbounded linear operators. This comes from the fact that when we consider the 
applications to partial dififerential equations then the operator associated with such 
equations are dififerential operators which are unbounded operators. 

There are many types of solutions of a dififerential or an integrodifferential equation 
in an abstract space. A form of solution which assumes the existence of derivatives 
of the solution up to the order of the equation is called classical solution. If we 
relax the condition of the existence of the derivatives upto lesser order, we get other 
type of weaker form of solutions. 

For establishing the existence and uniqueness of strong solutions in the case of short 
memory we use the method of semidiscritization. Using this method we may not 
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only establish the existence, uniqueness and regularity of solutions but obtain the 
approximate solution also in the process. This method so far has been successfully 
applied to various classes of linear as well as nonlinear parabolic and hyperbolic prob- 
lems. We replace the time derivatives by their corresponding difference quotients 
in the time dependent problems. This gives rise to a system of time-independent 
equations. These systems are guaranteed to have unique solutions, giving us the ap- 
proximate solutions, by the theory of accretive operators. After establishing some 
a priori estimates for the approximate solutions, one proves the convergence of the 
approximate solutions to the unique solution of the problem under consideration. 
The semigroup of operators are obtained as solutions of initial value problems for a 
first order equation in a Banach space. Most of the theory deals with a single first 
order equation. The reason for this is that higher order equations can be reduced to 
first order system and then by changing the underlying Banach space one obtains a 
single first order equation. 

Our aim in the present work is to consider the various physical problems of vis- 
coelastic solid models with short and long memory. As in the case of short memory 
of viscoelastic materials the problem is reformulated as a second order equation in 
a Banach space and the method of semi-discretization in time is applied to estab- 
lish the existence and uniqueness of strong solutions. In case of long memory the 
problem is reformulated as a single abstract Cauchy equation in a product Hilbert 
space. In the abstract Cauchy equation the operator is proved to be the infinitesi- 
mal generator of contraction semigroup in the product Hilbert space, and then the 
theory of semigroups is used to obtain the existence, uniqueness and asymptotic 
behavior of solutions of the problem. We also consider the exponential stability of 
the semigroup generated by the operator appearing in the abstract Cauchy problem. 
After obtaining the solutions for viscoelastic system of long and short memory, the 
regularity of such solutions is established in a Banach space. We also analyzed the 
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Fadeo-Galerkin approximation of solutions to the integrodifferential equations. We 
first consider an approximate equations and the existence of a unique solution to 
this approximate equation. After proving some estimates for the solution of the 
approximate equation, we establish the convergence of the approximate equation 
and its solution to the equation under consideration. 
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Chapter 1 


Literature Survey and Outline of 
Thesis 

1.1 Introduction 

The work presented in this thesis deals with the wave propagation and heat con- 
duction in the materials with memory effects. The mechanical behavior of solids 
is described by various models in continuum mechanics. In the present work we 
mainly concentrate on viscoelastic models. In Elasticity, an instantaneous response 
of material is treated as an adequate approximation. Finite elasticity is a part of 
elasticity theory and does not put any a priori limitations on displacements and 
strains. Synthetic rubber, rubber-like polymers and elastomer, biological tissues, 
metals and alloys under high pressure are some of the examples of elastic media 
with large deformations. 

Viscoelasticity may be considered as a theory taking into account of the effect of 
the entire deformation history on stresses at current instant. Finite viscoelastic- 
ity deals with nonlinear behavior at large deformations. Rubber-like polymers and 
plastics, metals at elevated temperatures, soils, road construction materials, biolog- 
ical tissues, food-stuflfe, polymeric melts and solutions are typical examples of such 
materials. These materials exhibit two different kinds of behaviors; liquid-like and 
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solid-like. Here, we are mainly concerned with conservation of momentum and en- 
ergy in viscoelastic media with finite deformations. 

Viscoelastic behavior is typical of a number of materials which are extremely impor- 
tant for applications. We refer to, for instance, Aklonis, MacKnight and Shen [1] for 
polymers and plastics, Alberola, Lesueur, Granier and Joanicot [2] for composites, 
Skrzypek [3] for metals and alloys at elevated temperature, Maccarrone and Tiu 
[4] for road construction materials, Deligianni,Marris and Missirlis [5] for biological 
tissues, Robert and Sherman [6] for food-stuffs. 

1.2 Theory of Viscoelasticity: A brief Introduc- 
tion 

Consider a rod-shaped specimen in its natural stress free state. Suppose we apply 
a tensile forces P to its ends at instant t = 0 and provide the stress (Jq = P/S 
where S is the cross sectional area. Immediately after application of tensile forces 
an instant strain eo arises in the specimen. In the case of pure elastic material, strain 
e{t) at instant t > 0 coincides with the initial strain cq- For an inelastic material 

p p 


Figure 1.1: A Specimen Under Uniaxial Load 

an additional strain e(t) > eo can be observed which increases monotonically with 
time, i.e., e(t) > 0, t > 0, where dot denotes the differentiation with respect to time. 

Definition 1.1 This phenomenon of growth of strain under a load uniform in time 
is called creep. 

The creep process is characterized either by the creep strain e^, given by 

edt) = e{t)-€o, 


( 1 . 2 . 1 ) 
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or by the creep compliance J which is the ratio of creep strain to the applied stress 
(7o, given by 

Jit) = ^ = iMufo. (1.2.2) 

<^0 CTo 

The growth of strain in time is a characteristic feature of both inelastic liquids as 
well as inelastic solids. For solid media, eit) — )■ Coo and e{t) — >• 0 as t oo with 
eoo < oo. For liquid media, eit) — >■ Coo as t — >• oo with €oo < oo corresponding to a 
Newtonian flow. The Newtonian liquid under stress cto uniform in time, the strain 
is proportional to time 

<«(«) = (1.2.3) 

where rj is called Newtonian viscosity and en(t) is called Newtonian strain. Taking 
into account (1.2.1), we may express 

€(t) = Co + + Cn(t)j (1.2.4) 

where cq is instantaneous elastic strain and edt) = edt) — edt) is the creep proper 
strain. Equation (1.2.3) is valid for both solid as well as for liquid media with rj = oo 
and 77 < oo, respectively. From (1.2.4), it follows that the creep proper strain tends 
to a finite limit for solids as well as for liquids. The division of media into solids 
and liquids is rather conditional as it depends essentially on temperature. Metals at 
room temperature are elastic solids whereas at elevated temperatures they demon- 
strate the creep typical of liquids. The strain €0 characterizes an instant response 
of material. In applications, media are classified accordingly with instant response, 
i.e., €0 7^ 0 and without instant response, i.e., €0 = 0 (cf. Giesekus [7]). 

Consider now another t3q)e of loading in which at instant i = 0 the specimen is 
stretched up to a length h and this new length remains fixed for all t > 0. Then 
the strain eit) = Cq for all t>0. For a purely elastic material the stress o-q arises in 
the specimen at < = 0 and remains unchanged during the loading. For an inelastic 
material the stress cr(t) = cto at i = 0 and ait) decreases monotonically as t increases 
and tends to a limiting 'ra.lue ergo as i -4 00. 
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Definition 1.2 The above phenomenon of decreasing of stress in time for a fixed 
strain is known as stress relaxation. 

If (7oo > 0 ) it corresponds to a solid and if a^o = 0 , it corresponds to a liquid. A 
new shape is created by stretching the specimen. A medium demonstrates the liquid 
type behavior if no additional stress arises, when a new shape without changes in it’s 
volume is provided. For uniaxial loading, the condition aoo == 0 means that asymp- 
totically ( for large time) the medium accepts its new shape without resistance, i.e. 
behaves as a liquid. 

The stress relaxation modulus E{t) is defined by 

E(t) = (1.2.5) 

Co 

For a linear elastic material, its elastic modulus E = a/e(j is inverse of the compliance 
J = e/cTo. In general E{t) is not the inverse of J{t). With the notion of creep and 
relaxation phenomena, we may define a viscoelastic material as follows. 

Definition 1.3 A material is called viscoelastic, if it exhibits both creep and 
relaxation phenomena. 

Again, this definition is rather conditional. Sometimes we refer to viscoelasticity 
when only one of these phenomena is taken into account. For example, some ap- 
plied problems of interest are concerned with either creep or relaxation. Creep and 
relaxation are processes which correspond to loads uniform in time when either 
stress or strain remains constant. In applied problems we encounter time-dependent 
loads. Under these loads the material exhibits a behavior which is not relevant to 
pure creep or relaxation. 


1.3 The Recovery Phenomenon 

When any of stress or strain is a step function 


a{t) = (JoUit), eit) = €on{t) 


(1.3.1) 
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then the creep and relaxation tests are the simplest experiments, where T-Lif) is the 
Heaviside function 


Hit) 


1 : t > 0 

0 : t<0 


(1.3.2) 


The above loadings cannot reveal all the characteristic features of the viscoelastic 
behavior, and more sophisticated programs are necessary. Several tests with time- 
dependent forces are employed (cf Szczepinski [8]), but the most convenient for 
applications are experiments with piece-wise constant loads. 

The recovery phenomenon is analyzed by using the following form of loading 


u{t) = 



t < 0 
Q<t<T 
T <t 


(1.3.3) 


The function a{t) and e(t) corresponding to equation (1.3.3) are plotted as follows. 
The creep recovery is the difference between the creep strain under loading and 
the creep strain when external forces and removed. We have two different ways for 
defining the creep recovery. 


Definition 1.4 According to Hadley and Ward [9], Ward and Wolfe [10], and Ward 
[11], the creep recovery equals the difference between strains in two specimen: one 
of them remains under the action of longitudinal forces and other is unloading : €r 
see Fig 1.2 (B) 


Definition 1.5 According to Schwartzl [12] and Gramespacher and Meissner [13], 
the creep recovery equals the difference between the strain at the current moment 
of time : see Fig 1.2(c) 


The second definition is very common for recovery strain. The typical dependencies 
of the total strain € = on time t for creep and recovery give rises to three 

cases: 

(i) for small stress, creep and recovery curves coincide; 

(ii) for large stresses, the “instantaneous” on short-time recovery is always greater 
than then “instantaneous” creep. 
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0 T 

(A) 




Figure 1.2: A Recovery Program For A Viscoelastic Specimen 


(iii) for any level of stresses, the creep and recovery strains increased in time mono- 
tonically. 


1.4 Nonlinearity of Material Response 

A number of elastic materials are rather hard and britle, and they demonstrate the 
linear response practically in the entire interval of strains before failure. Rubbers 
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Figure 1.3: Basic Rheological Elements 


and biological tissues are the only exceptions which permit large deformations up 
to several hundred percent. 

In viscoelastic media only a few materials demonstrate the linear response, whereas 
the majority of viscoelastic media are physically nonlinear. 

For linear materials, the ratio €c(t)/cro which determines the creep compliance J{t) 
should be independent of the applied stress ag. 

1.5 Differential Models with Small Strains 

Viscoelastic materials have both elastic and viscous features. In rheology, elastic 
and viscous properties are described by springs and dash-pots respectively. Fig 1.3 
gives the schematic representation of the basic rheological elements. 
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A linear elastic spring follows the Hooks law 

<7 = Ee, ( 1 . 5 . 1 ) 

where E is called Young’s modulus. For a nonlinear elastic string, we have 

a = $(e), (1.5.2) 

where $ is a smooth function. Viscous properties are modeled by dash-pots. In this 
case, the linear dependence of stress on the rate of strain is given by the Newton 
law. 


cr = T}e, (1.5.3) 

where 77 is called Newtonian viscosity. Again, in the nonlinear case, we have 

cr = ^r(g) (1.5.4) 

with ^ a smooth function. For basic rheological models, its natural to combine 
springs and dash-pots in a group. Combining in series leads to Maxwell model 
and the parallel combination gives rise to Kelvin- Voigt model. For more general 
rheological models we combine the basic elements (springs and dash-pots) as well 
as the Maxwell and Kelvin- Voigt elements in parallel and in series. 

If we take a Maxwell model with an elastic spring in parallel, we obtain the standard 
viscoelastic solid known as Zener model. The standard viscoelastic solid is one of 
the basic models for the theoretical analysis as this model describes both creep and 
relaxation processes and it can be represented both in the differential and integral 
forms. 

We now derive differential equation for a standard viscoelastic solid, for this we use 
the following rules. 

(i) If the elements are connected in parallel, their strains coincide and the total stress 
is the sum of stresses in separate elements. 

(ii) If the elements axe connected in series their stresses coincide and the total strain 
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Figure 1.4: The Maxwell Model 


equals the sum of strains in separate elements. 

Applying these rules to Maxwell model, we have 

e = €e "h €v (1.5.5) 


where Cg and are the strains in the spring and in the dash-pot, respectively. 
Differentiating (1.5.5), we get 

€ = eg + €«. (1.5.6) 


Assuming the linearity of the elements and using (1.5.1) and (1.5.3), we obtain 

(1.5.7) 

For Kelvin- Voigt model presented in Fig. 1.3, we have 

a = CTg-l-crt,, (1.5.8) 


a a 

€ = — -I — . 
E T] 


where Cg and <j„ are the stresses in spring and dash-pots, respectively. Assuming 
the linearity of material response, we get 


cr = Ee + rje. 


(1.5.9) 
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Figure 1.5: The Kelvin- Voigt Model 


Now we derive the constitutive equation for standard viscoelastic solid presented in 
Fig. 1.6. We denote by E and rj the Young modulus and the Newtonian viscosity 
for the Maxwell element and Ei the Young modulus of the additional spring. Thus, 
we have 

(7 = <7e -j- Cjn, (1.5.10) 

where ag and am are the stresses in the spring and in the maxwell element, respec- 
tively. Using (1.5.1) and (1.5.10), we get 

(Tm = a-a^ = a~Eie. (1.5.11) 

FVom (1.5.11) and (1.5.9), we get the constitutive equation for the viscoelastic solid 
as 

i + - = (l + §)£ + ^e. (1.5.12) 

E r) E rj 

For a more general rheological model with an arbitrary number of springs and dash- 
pots, the corresponding constitutive equation takes the form 

da . dP'a , t} , _u p 


(1.5.13) 
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where m and n are positive integers and Ak and Bi are material parameters. For 
m = n, Eq. (1.5.13) may be written in integral form. We prove this for the case of 
standard viscoelastic solid, i.e., m = n = 1. In this case we have 


Aid' + Aq(7 — Bic + Bq6, 


(1.5.14) 


If Ai = 0, the system is said to have short memory. If Ai 7^ 0 the system is said have 
long memory. In the case of long memory, the constitutive equation (1.5.14) may be 
rewritten in an integral form. We may rewrite it as first order linear nonhomogeneous 
equation in a, 


+ (1.5.15) 

Using the fact that <7(0) = 0 and that the linear first order nonhomogeneous equation 


^(t) + P{t)y(t) = Q{t), 


(1.5.16) 


has the solution 

2/W 




l/(^o)+ [ 

JtQ 


P(u)du 


Q{s) ds 


(1.5.17) 
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where to is any arbitrary constant, we have 

a{t) = C [Bieis) + Boe(s)] ds. (1.5.18) 

■^1 Jq 

Using integration by parts and the fact that €(0) = 0, we get 

f ds = /” ds. (1.5.19) 

Jo Al Jq 

Using (1.5.19) in (1.5.18), we obtain 

<^it) = f^€(t) + |^(l - ^ e(s)e-(^“/^^)(*-*) ds. (1.5.20) 

Now, the one dimensional equation of motion in the presence of only external force 
/ is given by 

putt = (^x + f (1.5.21) 

where u is the displacement and p > 0 is the material density. In the case of short 
memory, i.e., = 0 in the constitutive equation (1.5.14), using the fact that e = Ux, 

we obtain 


P^tt — CJqUxx “b ^x'^xxt "b f 


(1.5.22) 


and in the case of long memory we obtain 

puuit,x) = Cuxx + D / R{t- s)uxx{s,x) ds + f{t,x) (1.5.23) 

Jo 

for some constants C and D. The initial and boundary conditions may be written 
as 


u(0,a:) = Uo{x), —ix,0) = Ux{x), xe{0,l), (1.5.24) 

u(t,0) = u{t,l) = m^ t€{0,T]. (1.5.25) 


Note that the general Dirichlet boundary conditions (1.5.25) can be reduced to 
homogeneous Dirichlet boundary conditions, i.e., u{t,0) = u{t,l) = 0 for t € {0,T], 
by changing u{t,x) to u{t,x) — {x/l)'^(t) — {{I — x)/l)<l>{t). 
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1.6 Creep and Relaxation Kernels 

We discuss basic properties of creep and relaxation kernels. Let us consider a vis- 
coelastic specimen in its natural (stress-free) state in the form of a rectilinear rod. At 
instant r > 0, unit tensile forces are applied to the ends of the rod. The longitudinal 
strain p(t, r) at moment t > r is presented as the sum 

p{t,T) = pQ{r)+ pi{t,T) (1.6.1) 

where po = is the instant strain caused^ and ei = p(i,r) is an additional 

(creep) strain caused by the material, E{r) is called current elastic modulus, 
and C{t, r) is called creep measure. Function C{t, r) is assumed to be sufficiently 
smooth and to satisfy the condition 

C{t,t) = 0 (1.6.2) 

Now at the instant t = 0 a time varying longitudinal load is applied to the specimen 
and let a{t) denote the longitudinal stress. Function cr(t) is continuously differen- 


tiable and satisfies the equality 


O 

11 

b 

(1.6.3) 

Definition 1.6 The function 


K(t,T) = S(t)^[^^^j+C(t,T)] 

(1.6.4) 


is called the creep kernel. 


Now using the creep kernel the constitutive equation for linear viscoelastic media, 
where as the strain €{t) at instant t caused by the stress history {o-(r) : 0 < r < i} 
equals the sum of the strains caused by the elementary strams may be written as 

^ A:(i,r)£r(r) dr] (1.6.5) 

where as the first term on right hand side determines the instant elastic deformation, 
while the other term determines the creep deformation. 

For a number of materials, the effect of aging is weak and it may be neglected. But 
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from mathematical point of view, it means the current elastic modulus E can be 
taken as a constant, E{t) — E, while the creep kernel K may be considered as a 
function of the difference t — r only, K{t,T) = K{t — T). For a non-aging viscoelastic 
medium, constitutive equation (1.6.5) written as 

^ cr{t) + K(t- r)cr(r) dr] 

= ^ crW + ^ Co{t - r)cT(r) dr] (1.6.6) 

where 

K{t) = Coit), Co{t) = EC{t) 

Equation (1.6.5) and (1.6.6) determine the strain e as a function of the stress cr 
are called Creep equations. By solving these equations. Relaxation equations are 
obtained which express the stress cr as a function of the strain e. For the aging 
viscoelastic medium we get 

a{t) = E{t)€{t) - f dT (1.6.7) 

The function x(i> '^) is presented in the form 

x(t,T) = E{r)^Q{t,T) (1.6.8) 

where the Q{t, r) is called relaxation measure, and satisfy the condition 

Q{t,t) = 0 (1.6.9) 

Definition 1.7 The function 

is called the relaxation kernel. 


Now Equation (1.6.7) may be rewritten as 


r(i) = E{t) ^e{t) - Rit, r)€(r) dr] 
r(t) = E[€(t)-^*i?(t-T)€(r)dr] 
= E [€(t) - ^ Qoit - r)e{T) dr] 


( 1 . 6 . 11 ) 

( 1 . 6 . 12 ) 
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from mathematical point of view, it means the current elastic modulus E can be 
taken as a constant, E{t) = E, while the creep kernel K may be considered as a 
function of the difference t — r only, K{t,T) = K(t—r). For a non-aging viscoelastic 
medium, constitutive equation (1.6.5) written as 

~ h - r)a{T) dr] 

where 

K{t) = Co{t), Coit) = EC{t) 

Equation (1.6.5) and (1.6.6) determine the strain e as a function of the stress a 
are called Creep equations. By solving these equations. Relaxation equations are 
obtained which express the stress cr as a function of the strain e. For the aging 
viscoelastic medium we get 

< 7 ( 4 ) = B(4)€(4) - j‘ dT (1.6.7) 

The function x(^) is presented in the form 

X(t,r) = E{T) + Q(t,T) (1.6.8) 

where the Q(t, r) is called relaxation measure, and satisfy the condition 

Q(T,r) = 0 (1.6.9) 

Definition 1.7 The function 

is called the relaxation kernel. 


Now Equation (1.6.7) may be rewritten as 


r(i) = E(i)[€(t) i?(t,r)e(r) dr] 
r(t) = E[e(i)-^ i?(t-r)e(r)dr] 
= E [e(t) - ^ Qo{t - T)€(r) dr] 


( 1 . 6 . 11 ) 


( 1 . 6 . 12 ) 
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where 

Rit) = -Qo{t), Qoit) = ^ (1.6.13) 

Definition 1.8 A measure Q{t, s) is called regular if it is twice continuously differ- 
entiable. If a measure Q{t, s) is only differentiable, and it’s derivatives, the relaxation 
kernel R{t,s), has an integrable singularity at t = s, the Q{t,s) is called weakly 
singular. 

1.7 Regular and Weakly Singular Kernels 

The simplest of regular relaxation measures corresponds to the standard viscoelastic 
solid and is given as 

Qo{t) = -x[l - exp{-^)^ (1.7.1) 

where x is the material viscosity, and T is the characteristic time of relaxation. It is 
assumed that 0 < x < 1 s,nd T > 0. Using equation ( 1.6.3) and the differentiating 
equation( 1.6.4) we have 

R(t) = |exp( - i) (1.7.2) 

where Ji(t) is the relaxation kernel. 

equation ( 1.6.4) has important features for applications: 

(i) Its simple mechanical interpretation is that, it can be modeled by the system 
consisting of two elastic springs and dash-post. 

(ii) It allows the creep kernel to be found in explicit form. 

(iii) Describes qualitatively the material response observed in experiments for both 
creep and relaxation. In particular, it implies finite positive limits for the creep 
and relaxation kernels. The power-law relaxation measure, Findly et al. (1989) and 
Rabotnov (1969) 

Qo(t) = -(|:)“ (1.7.3) 

where a G (0, 1) and T > 0 are material parameters. Using equation (1.6.13) and 
(1.7.3), we have 

al 



= vJa-i{t) 


(1.7.4) 
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where 77 = Q;(a)!t““, Ja{t) is the Abel Kernel. 

In the next chapter, we shall consider the abstract formulations of (1.5.22) and (1.5.23) 
in a Banach space. These abstract formulations allow us to concentrate on the main 
features of the problem which are common to a class of such problems. Thus, any 
result proved for a prototype of problem in a class can be used for all the problems 
belonging to the class. 


1.8 Literature Survey 

Our aim in the present work is to consider the abstract formulations of various 
classes of problems to which (1.5.22), (1.5.23) and related problems are members. 
By considering their abstract formulations we may concentrate more on the main 
properties which remain invariant for the whole class. Also, it gives us clear idea of 
the basic features of a problem and helps in the useful generalization. The rich the- 
ory of functional analysis may be used to establish wellposedness and convergence of 
various type of solutions to the problems under consideration. Our main tools from 
the functional analysis are the semigroup of operators in a Banach space, theory of 
monotone (accretive) operators in a Banach space, the method of semi-discretization 
in time. 

Consider, for instance, the differential equation 

^ = Au{t), i>0 (1.8.1) 

at 

u(0) = Wo 

where u : [0, 00) ->■ R, and A is a real constant. We know that the solution of (1.8.1) 
is given by 

u{t) = t>0. 

If we define the maps T{t) : R — > R given by 

T{t)r = re^\ t > 0, 

then T{t), i > 0 is a family of bounded linear operators from R into R satisfying 
the following. 
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(a) Since T(0)r = r, T(0) = /, the identity map. 

(b) Since 


T{t + s)r 


re^t+s) ^ 


= (re^®)e^* = T(t)(re^*) 
= r(t)T(s)r, 


( 1 . 8 . 2 ) 


we have T{t + s) = T(t)T(s) for all nonnegative t and s. 

If in (1.8.1) we have u : [0, oo) — f R”, Uq € R” and A is an n x n matrix, then also 
the solution has the form 

u{t) = e*'^UQ, 


where 


e 


tA 



k\ ’ 


which is majorized by the exponential series 


OO 


E 


<‘ 11 ^ 11 * 

it! 


(1.8.3) 


More generally, if X is a Banach space and u : 
a bounded linear operator from X onto X, the 
problem. 


du 

dt 


(t) = Au{t), u(0) 


[0,oo) -y X, uo € X and A is 
solution of the abstract Cauchy 


= uo, 


is given by 


u{t) = 


where is defined similarly as in (1.8.3). Again, in the general situation, if we 
define for t > 0, T(t)x = e*^a: for every x E X, then u{t) = T{t)uo satisfies the 
abstract Cauchy problem. The family {T{t), t > 0} of bounded linear operators in 
X has the properties (a) and (b). Such a family is called semigroup of botmded 
linear operators in X. The operator A is called the infinitesimal generator of 
the semigroup. Thus if we know the semigroup generated by A then we know the 
solution of the Cauchy problem. The main importance of the theory of semigroup 
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appears when we consider A as unbounded linear operator. This comes form the 
fact that when we consider the applications to partial differential equations then the 
operator A is associated with a differential operator. 

The theory of semigroup is dealt with in many books on functional analysis. One 
of the books with extensive work on the subject is the classical book by Hille and 
Phillips [14]. We also refer to other books of Butzer and Berns [15] Dunford and 
Schwartz [16] Ladas and Lakshmikantham [17], Martin [18], Yosida [19] and Gold- 
stein [20]. After the result of Hille- Yosida theorem in 1948, which gives the char- 
acterization of the infinitesimal generator of a contraction semigroup, there was a 
rapid development in the theory of semigroup. By now, it is an extensive math- 
ematical subject with substantial applications to many fields of analysis. Most of 
the classical results for semigroup of bounded linear operators in a Banach space 
were generalized to equi-continuous semigroup of class Co in locally convex spaces 
by Yosida [19]. The strongly continuous semigroup are considered in Pazy [21]. For 
the theory of nonlinear semigroup in Hilbert and Banach spaces, we refer to Barbu 
[22], Yosida [19]. 

In the third chapter we consider the case of short memory. We use the method 
of semi-discretization in time to establish the existence and uniqueness of strong so- 
lutions. The method of discretization in time was introduced by Rothe [23] and 
developed by a group working at Prague, see Rektorys [24], Necas [25], Kacur 
[26], Kartsatos and Zigler [27]. See also, Bahuguna and Raghavendra [28, 29] and 
Bahuguna, Pani and Raghavendra [30] and references cited therein. This method 
consists in dividing the time axis and replacing the time derivatives by their cor- 
responding difference quotients in the time dependent problems. This gives rise 
to a system of time-independent equations. These systems are guaranteed to have 
unique solutions, giving us the approximate solutions, by the theory of accretive 
operators. After establishing some a priori estimates for the approximate solutions, 
one proves the convergence of the approximate solutions to the unique solution of 
the problem under consideration. As mentioned above, semigroup of operators are 
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obtained as solutions of initial value problems for a first order equation in a Banach 
space. Most of the theory deals with a single first order equation. The reason for 
this is that higher order equations can be reduced to first order systems and then 
by changing the underlying Banach space one obtains a single first order equation. 
This idea we have used in the fourth chapter. The initial works on linear viscoelastic 
materials with short as well as long memory is considered by Glowinski, Lions and 
Tremolieres [31], Duvaut and Lions [32] and others. The further study has been 
considered by Dafermos [33, 34], Day [35], Desch and Miller [36] - [37], Hannsgen 
and Wheeler [38, 39], Fabrizo and Lazzari [40], Liu and Zheng [41, 42] and others. 
We consider a more general problem in the fourth chapter. In the fifth chapter, we 
consider the regularity of the solutions whose existence and uniqueness is consid- 
ered in the earlier chapters. In the last two chapters we consider the problem of 
heat conduction in material with memory and use the idea of fractional powers of a 
dissipative operators which was developed by Komatsu [43] Kato [44]. This allows 
us to consider more general nonlinear forcing terms and establish the existence and 
uniqueness of classical solutions. 


1.9 Outline of Thesis 

In this section we briefly mention the topics covered in the chapters two to seven of 
the thesis. 

Chapter Two; In this chapter we mention the basic notions, definitions and results 
from the literature. We first give the idea of abstract formulations of various prob- 
lems and their importance. We then consider the theory of semigroup of bounded 
linear operators in a Banach space. We then briefly describe the method of semi- 
discretization in time. 

Chapter Three: This chapter deals with the case of short memory of viscoelastic 
materials. We reformulate the problem as a second order equation in a Banach space 
and apply the method of semi-discretization in time to establish the existence and 



1.9 Outline of Thesis 


20 


uniqueness of strong solutions. We also make some remarks about the extension of 
results for the existence and uniqueness of classical solutions 

Chapter Four: This chapter considers the case of long memory. We reformu- 
late the problem as a single abstract Cauchy equation in a product Hilbert space. 
We prove that the operator appearing in the abstract Cauchy equation is the in- 
finitesimal generator of a contraction semigroup in the product Hilbert space. We 
then use the theory of semigroup to obtain the existence, uniqueness and asymptotic 
behavior of solutions of the problem. We also consider the exponential stability of 
the semigroup generated the operator appearing in the abstract Cauchy problem. 

Chapter Five: In this chapter we establish the regularity of the solutions for 
an abstract integrodifferential equation in a Banach space. The problem considered 
in this chapter includes the problems considered in the earlier chapters. Thus, the 
results established in this chapter prove the regularity of solutions obtained in the 
earlier chapters. 

Chapter Six: In this chapter we consider the Fadeo-Galerkin approximation of 
solutions to the integrodifferential equations analyzed in Chapter Six. We first 
consider an approximate equation and the existence of a unique solution to this 
approximate equation. After proving some estimates for the solution of the approx- 
imate equation, we establish the convergence of the approximate equation and its 
solution to the equation under consideration and its solution, respectively. 


The relevant references are appended at the end. 




Chapter 2 

Preliminaries and Basic Results 


2.1 Introduction 

In this section we mention some of the results required in the ensuing chapters. 
Most of the results are stated here without proofs but with proper references. 


We first consider the idea of an abstract formulation and its importance. We then 
consider various tools to study such abstract formulations. Some of these tools are 
the theory of semigroup of bounded linear operators in a Banach space X, mono- 
tone (or accretive) operators, method of semi-discretization in time. The notion 
of a semigroup iS(t), t > 0 of bounded linear operators from X into X and its in- 
finitesimal generator A are closely related with various partial differential equations. 
The generator A is the partial differential operator appearing in the partial differ- 
ential equations which model many physical phenomena. A typical example of this 
operator is 


A = 


-A = 


-E 


dxj 


which is n-dimensional Laplacian and appears in the heat conduction or reaction- 
diffusion equations and wave equations. 



2.2 Abstract Formulation 


22 


2.2 Abstract Formulation 

In order to explain the basic idea, we consider the following initial boundary value 
problem for the heat equation: 

~ f2x(0,T), (2.2.1) 

i=l 

M = 0 on d^x [0, T] 
u(rc, 0) = t4o(a:) on 

Here Q is a bounded domain in R" with sufficiently smooth boundary dO.. The 
functions / and uq are given and u is to be found. 

The generalized problem corresponding to (2.2.1) is obtained by multiplying (2.2.1) 
by a function v in the space C^(n) of infinitely differentiable functions having com- 
pact support in and integrating by parts. Here by support of a function cj) £ C{Q), 
we mean the closed set 

K = {x € f2 1 4>{x) ^ 0} 

We thus obtain 

u{x,t)v{x) dx 4- I^Diu{x,t)Div{x) dx 

i=i 

= f f{x, t)v{x) dx, X G n (2.2.2) 

Jn 

u{x,0) = ■uo(a;)5 xeQ,, 
where x = {xi,X 2 ,- • • , Xn) and A = dj dxi. We consider 

V = Wo^’^(^^), ff = L^(n) 

where W”‘’P(n) and W2”’*’(n) are the usual Sobolev spaces. (Cf. Pazy [21] for 
definitions and other details). We observe the following main points, 

(a) the treatments to the space variables x = (xi, X 2 , • • • , x,i) and the time variable 
t are different, (b) For fixed time t, the function 



X I ^ ti(x, t) 
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is an element of the Sobolev space V. We briefly denote this element by u{t), i.e., 
for each t € [0, T] 

u(t) € V. 

Thus, if we vary t in [0, T] then we obtain a function 

1 1 — ^ u(t') . 

Hence, from the real function {x,t) i — u{x,t), we obtain a Banach space V valued 
function 1 1 — > u(t). 


We now write the equation (2.2.2) as 


d 


+ a{u{t),v) 

- on (0,r), for all v€V, 

u(0) = Wo G H. 


(2.2.3) 

(2.2.4) 


We have to find the function 1 1 — > u(t) € V for all t € [0,T]. For all w, v, we define 

a(w,v) = / '^diw{x)Div{x) dx, 

i=i 

(/(*)> ■y) = [ f{x,t)v{x)dx. 

Jn 

Thus we see the natural need of considering two spaces H and V. H is required for 
the time derivative and V is needed for the space derivatives in the Laplacian 


-A 


V — 
^dxY 

*=i * 


In the present case, we have the embedding V ^ H \b continuous and V is dense 
in H. This leads us to the evolution triplet 


VCHCV*. 

The time derivative d/dt is considered as the generalized derivative on (0,r). For 
this reason equation (2.2,3) is required to be satisfied only almost everywhere (writ- 
ten in short as a.e.) on (0,T). For the function 1 1 — ¥ u{t) we choose the Sobolev 
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space 


W^'^{0,T-,V,H) = {ueLP{0,T-,V)\u' eL%0,T-,V*)} 


The corresponding norm is given by 

rT 




= [I IWWII 


X 1/2 

y dtj + 


/ .T ^ 1/2 

^ \\u'{t)ry.dt 


where u' is the generalized derivative of u. The generalized derivative of u is defined 
in such a way that 

d 

<u'{t),v>y = 3 T(w(i),w)H for all veV 
at 


To write the equation (2.2.3) as operator equation we introduce the operator A : 
V — > V* and the functional F{t) € V* through the equations 

< Aw^v >y = a{w,v) 

<F{t),V>y = {f{t),v)H 

for all u € y. We then obtain 

< u'{t) + Au{t) - F{t), V >y = 0 on (0, T), for all u € V 

which leads to the operator equation 

u'{t) + Au{t) = F{t), a.e.t€(0,r) (2.2.5) 

'ti(O) = uq ^ H. 

The equation (2.2.5) is an abstract first order equation known as first order evo- 
lution equation. One point we observe here is that u{t) € V and u'{t) G V* and 
we necessarily need to consider the evolution triplet V Q H Q V* in this abstract 
formulation. 


There is still another approach in whidi we need to consider only one space, namely 
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H = Z/^(Q) provided f{t) € H. For this we write the original problem (2.2.1) in the 
following form 

u'(t) + AHu{t) = f{t), for all t E {0,T) (2.2.6) 

u{0) = uo (2.2.7) 

where the operator Ag is the H space realization of A: 

DiAH) = {ueV\AueH}, 

Ahu = Au for u E D{Ah) 


2.3 Semigroups and their Generators 

Definition 2.1 A one parameter family S{t), t > 0 of bounded linear operators 
from X into X is called semigroup of bounded linear operators on X if 

(i) 5(0) = J where I is the identity operator on X, 

(ii) S{t + s) = S{t)S{s) for every t,s >0. 

The linear operator A whose domain is defined by 


D{A) 

f VI- S{t)x-x 1 

40 t J 

(2.3.1) 

and the action of A on x 

E D{A) is given by 



Ar = 

40 t 

(2.3.2) 


is called the infinitesimal generator of S{t). 
Example Let S{t) : R — ^ R be given by 


S{t) = e^^xo, t>0, xo G R. 

Then S{t) is a semigroup of bounded linear operators on R and its infinitesimal 
generator is A = XL | 

Exsimple Let S{t) : R” — > R", be such that 

S{t) = e*®Xo, t>0, Xq € R" 
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where B = [6ij] is an n x n matrix, bij are constants. Then S{t) is a semigroup of 
bounded linear operators on X R” and its generator is A = B. | 

Example Let S(t) : X — > X be given by 

S{t) = e^^xo, t > 0, xo € R 

where B : X — > X is a bounded linear operator. Then S{t) is a semigroup of 
bounded linear operators on X. and its generator is A = B. | 

Definition 2.2 The semigroup S{t) is called uniformly continuous if 

lim||5(t)-/|l = 0. (2.3.3) 

We note that the semigroup S (t) in each of the above examples is uniformly contin- 
uous. In fact, we have the following result. 

Theorem 2.1 A linear operator A is the infinitesimal generator of a uniformly 
continuous semigroup if and only if A is a bounded linear operator. Furthermore, 
the semigroup S{t) generated by a bounded linear operator A : X — y X is unique 
and 

(a) There exists oo >0 such that |l<S'(i)|| <e^. 

(b) The map ty-yS{t) is differentiable in norm and 

= AS{t) = S(t)A. 

For applications to partial differential equations, we can not expect the infinitesimal 
generator of a semigroup to be a bounded linear operator. Therefore, the class of 
uniformly continuous semigroups is of little use for this purpose. Thus, we need to 
consider other type of semigroups. 

Definition 2.3 A semigroup S(t), t > 0 of bounded linear operators on X is called 
strongly continuous or Co - semigroup if 

limS'(t)x = X 

40 



for all X € A. 
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For Co-semigroups, we have the following main result. 


Theorem 2.2 Let S(t) be a Co-semigroup and let A be its infinitesimal generator. 
Then 

(a) There exist constants M >1 and u >0 such that 

|15(t)I| < for t>0. 


(b) For X E X, 

lina - j S{s)x ds = S{t)x. 

(c) For X E X, 

J S(s)x ds E D(A) and A S{s)x 

(d) For X E D{A), S{t)x E D{A) and 


= S{t)x — x. 


dt 


S{t)x = AS{t)x 
= Sit)Ax. 


(e) For X E D{A), 

S{t)x - S{s)x 




dr. 


One of main results in the theory of semigroups of bounded linear operators is 
the Hille-Yosida theorem which was established in 1948 by Hille [45] and Yosida 
[19]. This result establishes a characterization of the infinitesimal generator of a 
semigroup S{t) with the property that 


115(t)l| < 1 for t>0. 


(2.3.4) 


Definition 2.4 If a semigroup t > 0 of bounded linear operators on X satis- 
fies (2.3.4), it is called Co-semigroup of contractions on X. 


For the characterization of the infinitesimal generator of a semigroup of contractions, 
we have the following result. 
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Theorem 2.3 (Hille-Yosida) A linear operator A is the infinitesimal generator of 
a Co-semigroup S{t) of contractions, t>0 if and only if 

(i) A is closed and D(A) is dense in X, 

(ii) the resolvent set p{A) of A contains [0, oo) and 

||{A/-.1)-‘|| < i 

for every A > 0. 

We note that if S{t) is a Co-semigroup with A its infinitesimal generator such that 
for some cj > 0, 


IlSWII < e"' (2.3.5) 

then Si{t) = e~‘^S{t) is a Co semigroup of contractions and its infinitesimal gener- 
ator is A — (vI. Now we consider the notion of dissipative operators. To define the 
dissipative operators we require duality map defined from X into the dual space X* 
of X, given by 

F{x) = {a;* & X* : <x*,x >= Hxlp = lla:*||^} (2.3.6) 

where < ., . > denotes the duality paring between the dual space X* of X and X. 
Clearly, duality map is a multi-valued or set valued map in general. In particular, 
when X is a Hilbert space then the duality map is single-valued. 

Definition 2^5 A linear operator A is called dissipative if for every x € D{A) 
there exists a x* such that 


Re < x*, Ax > < 0. 

Following theorem gives the diaracterization of dissipative operators. 
Theorem 2.4 A linear operator A is dissipative if and only if 

Alixil < 11(A/-A)x|l 


for all X € D{A) and A > 0. 
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Another useful characterization of the infinitesimal generator of a Co-semigroup of 
contraction is given by the following theorem due to Lumer and Phillips [46]. 

Theorem 2.5 (Lumer and Phillips) Let A be a linear operator with dense domain 
D{A) in X. Then we have the following. 

(a) if A is dissipative and there exists a Xq > 0 such that the range R{XoI — A) = X, 
then A is the infinitesimal generator of a Co-semigroup of contractions on X. 

(b) if A is the infinitesimal generator of a Co-semigroup of contractions on X then 
R{XI — A) = X for all X> 0 and A is dissipative. 

Moreover, for every x € D{A) and every x* € F(X), Re < x*,Ax > < 0. 


As pointed out earlier that what we know is the operator A associated with some 
partial differential operator and the problem is to find the semigroup generated 
by it. For a bounded linear operator B, we already know that the semigroup is 
uniformly continuous and given by e*^. The main task is to find the representation 
of the semigroup generated by a unbounded linear operator, if there is any. For this 
purpose we use the complex operator calculus. For a bounded linear operator B we 
have another representation given by the following theorem. 


Theorem 2.6 Let B be a bounded linear operator. If'f> |lBj|, then 

^ 74-100 




1 pj+too 

= - B)-'^dX. 

27ri J y—i(X3 


(2.3.7) 


The convergence in (2.3.7) is in the uniform operator topology and uniformly in t 
on bounded intervals. 


The following theorem gives the sufl&cient conditions for an operator A to be the 
infinitesimal generator of a Co-semigroup and a similar representation as in (2.3.7). 

Theorem 2.7 Let A be densely defined operator in X satisfying the following con- 
ditions. 

(i) For some 0 < 5 < 7r/2, p{A) D Sj = {A : jor^ A| < (7r/2) -H 5} U {0}. 

(a) There ansts a constant M such that 


M 
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for A G Ej with A > 0. 

Then, A is the infinitesimal generator of a Co-semigroup S{t) satisfying |jiS(t)|| < C 
for some positive constant C and 

S{t) = /" e^\XI - A)-'^dX 

ZTCl Jp 

where T is a smooth curve in T!,s running from ooe~'^ to ooe^ /or 7r/2 < 9 < 7 r/ 2 -\- 5 . 
The integral converges for t > 0 in the uniform operator topology. 

2.3.1 Semigroups of Compact Operators 

We recall that a linear operator T : X — y X is called compact operator if for 
any bounded subset D C X, T{D) is relatively compact in X, i.e., the closure T{D) 
of T(D) in X is compact. 

Definition 2.6 A Co-semigroup S{t) on X is called compact for t > to if for every 
t > to, S{t) is a compact operator on X. S(t) is called compact if S{t) is compact 
for t > 0. 

Note that we can not allow S{t) to be compact for t > 0 because in that case the 
identity operator will be compact and hence the space X has to be finite dimensional. 
Also, if S{to) is a compact operator on X for some to > 0? then S{t) = S(t—to)S{to) 
for t > to and therefore S{t) is compact as S{t - to) is bounded. Hence if S{t) is a 
compact operator for t > to if and only if 5(to) is a compact operator. We have the 
following main result for a compact semigroup. 

Theorem 2.8 Let S{t) be a Co-semigroup. If S{t) is compact fort > to, then S{t) 
is continuous in the uniform operator topology for t> to. 

The following result gives the characterization of the infinitesimal generator of a 
compact semigroup. 

Theorem 2.9 Let S{t) be a Co-semigroup and let A he its infinitesimal generator. 
S{t) is a compact semigroup if and only if S{t) is continuous in the uniform operator 
topology for t > 0 and {XI — A)~^ is compact for X € p{A). 
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2.3.2 Analytic Semigroups 

This class of semigroups on X extends the notion of a Co-semigroup S(t) on X 
where the parameter t lies on the nonnegative real axis to the semigroup S{z) where 
the parameter z lies in a sector 

A = {z : 4>i < arg 2 < 02} (2.3.8) 

which contains the nonnegative real axis and is analytic on A. To be precise, 
we have the following definition and theorem. 

Definition 2.7 Let A be as in (2.3.8) and for z € A S{z) be a bounded linear 
operator on X. The family S{z), z £ A is m analytic semigroup on X in A if 

(i) z S{z) is analytic in A, 

(ii) 5(0) = I and lim 2 _,.o, 2 gA 'S'( 2 ;)x = x for every x £ X, 

(iii) 5(^1 -f Z 2 ) = S{zi)S{z 2 ) for Zi,Z 2 £ A. A semigroup S{t) is called analytic if 
it is analytic in some sector A containing the nonnegative real axis. 

We note that the restriction of an analytic semigroup to the nonnegative real axis 
is a Co-semigroup. When can we extend a Co-semigroup to an analytic semigroup 
in some sector containing nonnegative real axis? Prom Theorem 2.7, we have that 
for a densely defined operator Ain X satisfying (i) and (ii) of Theorem 2.7 is the 
infinitesimal generator of uniformly bounded Co-semigroup S(t). Furthermore, in 
this case, S{t) can be extended to an analytic semigroup in the sector 

As = {z : |arg z\ < 5} 

and 1 1 5 ( 2 :) II is uniformly bounded in every closed subsector As>, S' < 5. We collect 
all these result in the following theorem. 

Theorem 2.10 Let T{t) be a uniformly bounded Co-semigroup. Let A be the in- 
finitesimal generator of S{t) and assume that 0 £ p{A). Then the following are 
equivalent. 

(a) S{t) can be extended to an analytic semigroup in a sector Aj and l|5(t)|| is uni- 
formly bounded in every closed subsector As> of As, S' <S, 


2.4 Method of Semi-discretization in Time 


32 


(b) There exists a constant C such that for every a > 0, r ^ 0, 

||((a-hir)/-4)~^|| < ^ 

Fl 

(c) There exist 0 < S < 7r/2 and M > 0 such that 

p{A) D S = |a : jar^ A1 < I -f 5| U {0} 

and 

||(A/-.4)-'|l < ^ 

for A E S, A ^ 0. 

(d) S{t) is differentiable for t > 0 and there exists a constant C such that 

||>ts(()|| < J, t>0. 

2.3.3 Fractional Powers of Operators 

In the case when —A is the infinitesimal generator of an analytic semigroup, we 
may define the fractional powers of A. This allows us to consider in the ensuing 
chapters the nonlinear maps involving differential operators also. As we know from 
Theorem 2.10 that if A is densely defined closed linear operator satisfying 

p(A)E+ = {A : 0 < |arg A| < tt} U V, (2.3.9) 

where V is a neighborhood of zero, and 

|I(A/ - .4)-‘|| < A £11+ (2.3.10) 

2.4 Method of Semi-discretization in Time 

In 1930 Rothe [23] suggested a method in which he converted a one dimensional 
parabolic initial-boundary value problem into a system of boundary value problems 
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for ordinary differential equations. He divided the time interval [0,T] into n subin- 
tervals j = 1,2, ...,n where tj = j.h and h = T/n, and replaced the 

partial derivative du/dt at (x,t) by the difference quotient 

u]{x) - 
h 

This idea was adopted in a large number of later works which treated more general 
parabolic and hyperbolic problems. For more details on the initial development 
of the method, we refer to the works of Rektorys [24]. For extension to nonlinear 
problems, we refer to Necas [47], Kacur [26], Kartsatos and Zigler [27] and references 
cited therein. 

In these works, estimates have been proved directly for the difference quotients. 
This requires global Lipschitz condition on the nonlinear forcing terms. Bahuguna 
et al. [48, 49, 50, 30] have modified the method so that one first proves the a priori 
estimates for the discrete points and then with the help of Local Lipschitz conditions 
one establishes the a priori estimates for the difference quotients. This allows us to 
consider more general nonlinear forcing terms. 



Chapter 3 


Viscoelastic Materials With Short 
Memory 


3.1 Introduction 


In this chapter we shall consider the case of short memory. We reformulate (1.5.22) 
as a second order abstract Cauchy problem in a Banach space. Our aim is to study 
the existence and uniqueness of strong solutions and our tool is the method of semi- 
discretization in time. This method has been successfully applied to various first and 
second order linear as well as nonlinear unsteady state problems. See, for instance, 
[28, 29, 30]. The method consists in dividing at each step n = 1, 2, • • •, the finite time 
axis [0,T] for arbitrary 0 < T < oo into subintervals — jh, h = T/n 

and replacing the time derivatives 


du 

dt 


5u] = 


u] - 


^ _ X 2 ,,n ^ U] -2u]_-^+u]_2 

at2 J A;2 


After establishing a priori estimates on these difference quotients one finally proves 
the convergence of the approximate solution 

lP{t) = + 
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to the solution of the problem. 


3.2 Abstract Formulation of the Problem 

In this section we describe the abstract formulation of the problem. We reformulate 
the problem in an abstract space so that our problem can be treated as one of a 
large class of such problems. Let X = 1/ (0,1), 1 <p < oo, and let u : [0,T\ X 
be given by 


u{t){x) = u{t,x), for all x E {0,1), t€[0, T]. 

Let L be the linear operator defined by 

D{L) = 

cPu , . 

Lu = — for uED{L), 
dx^ 

where and for an open subset of the n-dimensional Euclidean 

space R”, are Sobolev spaces (cf Pazy [21] for the definitions). Under reason- 
able assumptions on /, <j) and and homogeneous Dirichlet boundary conditions, 
Eq. (1.5.22) can equivalently be written in the operator form in X as 

+ aL^ + hLu(t) = /(t), (3.2.1) 

ii(0) = iJo, ^ = Vi, 

where / :‘[0,T] ->■ X given by 

f{t){x) = cf{t,x), for all a;G(0,Z), te[0,r] 

and a > 0, b and c are some constants. For notational convenience, we replace u 
and f hy u and /, respectively. Transferring other terms on the right in (3.2.1), we 
obtain 

t»(0) = Uo, ^ = Uu 


(3.2.2) 
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where A = Aq/ + aL for some suitable positive constant Aq and for some constants 
a and /3. Here I is the identity operator on X. Also, D{B) = D{A) with 

Bu = cxAu + C 2 U for u e D{A) 

for some constants Ci and C 2 . 


We assume the following properties of the linear operators A, B. 

(PI) The linear operator A is m-accretive, i.e., for every u € D{A) and A > 0 


A|lu|| '< ||Au + Ati|| 

and for every A > 0 

il(H-AA) = X, 

for A > 0 where R(.) is the range of an operator. 
(P2) (/ 4- is compact. 

(P3) The linear operator B satisfies 


l|Bwl| < C[ ||Au|l + |]«|| ^ 


for u G D{B), where C is a positive constant. 

If the operator —A is the generator of a Co-semigroup of contractions in X then 
it follows from Lummer-Phillips theorem (see 1.4.3 in [21]) that A satisfies (PI). 
The property (P2) follows from the compact embedding of Wo’^(0, i) in 1/(0, Z). 
Property (P3) is straightforward to check. 

By a strong solution to (3.2.2) we mean a function u : J — y X such that du/dt is 
absolutely continuous on J into X, u(0) = Uo, du{0)/dt = Ui and (3.2.2) is satisfied 
a.e. on J. 

Our aim is to prove the following main result for the initial data {Uq, Ui) in D{A) x 
D(A). 

Theorem 3.1 There exists a unique strong solution to (3.2.2) on J for every pair 
(Co, Cl) inD{A)xD{A). 
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3.3 Proof of Theorem 3.1 


In this section, we prove Theorem 3.1 with the help of Lemmas 3.3 to 3.6, stated 
and proved below. For a positive integer n, we consider the discretization 

t] = j.k, k = j = l,2,---,n; 
n 

of the interval J. We call u” an approximate solution and set 


Uq = Uo, 

= Uo-kUi, (3.3.1) 

= k^[f{0,Uo,Ui) + AUQ]-2kUi + Uo, 


for every n and discretize (3.2.2) as follows: 

+ A5u] + Bu]_, = /;, i = 1, 2, • • • , n; (3.3.2) 


where 


5u^ 


u] - 2u]_i + u ]_2 

A:2 

u] - u"_i 
k 


(3.3.3) 


+ + j = l,2,---,n. 

The existence and uniqueness of m” € D{A), j = 1, • • • ,n; n > 1, is a consequence 
of m-accretivity of A. For notational convenience, we suppress the superscript n 
and use C to denote a generic positive constant independent of the discretization 
parameters j, k and n, i.e., we allow C to have different values in the same discussion. 
We shall use the following lemma which is due to Sloan and Thomme [53] to get the 
estimates in Lemma 3.3. 


Lemma 3.2 Let {tOn} be a sequence of nonnegative real numbers satisfying 


n— 1 


Wji ^ OCji d" n ^ 0, 


h=0 


where {o!n} is a nondecreasing sequence of nonnegative real numbers and ^ 0- 
Then 

< OJnexp ( ), n>0 


w„ 


s.^i=0 
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The following lemma gives the uniform bounds for the difference quotients {S^Uj} 
which in turn imply the uniform bounds of {5wj}. This is essential for the conver- 
gence of the approximate solutions to the solution of (3.2.2). 

Lemma 3.3 For all n and j = 1, 2, • • • , n 

l|i\.|| < a 

Proof. We recall here' that by C we shall denote a generic positive constant inde- 
pendent of j, k and n. We allow C to have different values depending only on the 
constants T, ||.4/7o||) ll-6t/o||, 11/(0 , C/q) etc. Therefore, the constants 

CT, Ce^'^, etc., will again be replaced by C. 

For / = 1 in (3.3.2), we have 

5\i + kA5\i = A - BUq - AU^. 


The accretivity of A implies that 

IliVll < c. 

Now, for 2 < / < n in (3.3.2) we have 

5'^Uj -t- kAS^Uj = — kB5uj-i + fj — fj-i. 


Again, the accretivity of A implies 

PXII < ||dVill + *l|B'5%-ill + ll/i-/»-ill- 

Reiterating the above inequality and using (P3), we get 

i-i j-i 

\\SW\ < C + CkJ^\\SW\ + CkJ2U^^jl 

t=l i=l 

Now from (3.3.2), we have 

Ipfaill < c + ||d"tii||, 

and for 2 < i < n, we have 

i-l i -1 

i|A5Uil| < C+\\5W\+CkJ^\\S^Ur\\+CkJ2U^'^r\\. 


r=l 


r=l 


(3.3.4) 
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In view of the fact that jk <T hr 1 < j <n, we have the estimates 

l|A5uill < C 1+ max l|(5^tirl| + l|i4^Ur|l • 

r=0 

Lemma 3.2 implies that 

||A5uj|l < C 1 + max l|(5^Wr|| , 1 < z < n. (3.3.5) 

Using the estimates of (3.3.5) in (3.3.4) we get 

r 7-1 

max ||5^ztr|| < C l + kS^ max 115^^11 • 

L i=0 J 

Again using Lemma 3.2, we get the desired result. 

Remark 1 The estimates of Lemma 3.3 imply that for all n and y = 1, 2, • • • , n, < 
C. 

Definition 3.1 We introduce sequences oi polygonal functions {U”} and {U"} from 
J into X, given by 

U (t) = Uj—i “h (t tj—i)5Uj.) t € [tj— 1, tjj) J ~ 1» 2, • • • , 7T., 

(t) = 5uj-i + {t- t € [tj.i ,tj], J = 1, 2, • • • , n. 

Further we define the sequences of step functions {A"} and {T”}, from {—k, T] into 
^,by 

X"(t) = Uo for tG(-A:,0], X^{t) = uj for 

Y^-it) = Ui for te(-k,0], Y'^it) ^ 6uj for 
for j = 1, 2, • • • , n. 

Remark 2 The sequences {U"(t) - X”(t)} and {W‘(t) - F"'(i)} converge to zero 
uniformly on J as n — Y oo. Moreover, the functions {U^} and {U"} are Lipschitz 
continuous on J with uniform Lipschitz constants. 
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Lemma 3.4 There exist subsequences {[/”*} of {[/"} and {V'^} of {F”} and Lip- 
schitz continuous functions u and v such that — )■ u and V”* — )• v in C(J, X) 

(with supremum norm) as m — > oo. Furthermore, 

du (Fu dv 

— = u on J and = — a.e. on J. 
dt dt^ dt 

Proof. The uniform bounds on <5^u/s in the estimates (3.3.5) imply that Adu/s 
(and therefore Auj^s) are uniformly bounded for all n and j = 1,2, - ■ ■ ,n. We use 
(P2) to conclude that there exist subsequence and {T*”} of {X”} and {T"}, 

respectively, and functions u and v defined from J into X such that 

X'”(t) ^ u{t), and ^ u(t) 

as m — > oo. 

Now, the uniform Lipschitz continuity of {X™} and {y*"} imply that the families 
and {y*^} are continuous and the convergence of {X’"} and {y"*} to u and 
V imply that {X'”(t)} and {y”^(t)} are relatively compact for each t in J. The 
Ascoli-Arzela theorem implies that 

U^—^u and 


in C{J, X) as m — > oo. 

Since {U'^} and {V"^} are Lipschitz continuous with uniform Lipschitz constants, 
it follows that u and v are Lipschitz continuous on J. Now, for each x* G X*, we 
have 

/■* dU^(s) 

<U'^{t),x*> = / <^=^^^,x*> ds+<Uo,x*> 

Jo dt 

<V'^{s),x*> ds+<Uo,x*>. 

Passing to the limit m — > oo, we obtain 

<u{t),x*>= f <v(s),x*> ds+<Uo,x*>. 

Jo 

Therefore dufdt = v on J. Hence, 



(Fu 

dt^ 


dv r 

— a.e. on J. 
dt 
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Lemma 3.5 For t e J, u{t), du/dt e D{A), AX^(t) Au{t), AY^{t) 
A{du{t)/dt) uniformly on J as m — y oo and Au{t), A{du(t)/dt) are weakly contin- 
uous on J. (By we denote the weak convergence in X.) 

Proof. Since X'^{t) — y u{t) and Y^{t) — y du{t)/dt uniformly on J as m — y oo 
and {AX^(t)}, are uniformly bounded, we use Lemma 2.5 due to Kato 

[54] to assert that for t e J, u{t), du/dt e D{A) and AX'^ Au{t), AY^{t) 
A{du{t)/dt) uniformly on J as m — y oo. 

To prove the weak continuity of Au{t) and A{du{t)/dt) on J, let {ts} C J be such 
that ts — y t E J as s — y oo. Then u{ts) — y u{t) and du{ts)/dt — y du{t)/dt as 
s — y oo. Since Au{ts) and A{du{ta) J dt) are the limit of and 

as m — y oo, {Au(ts)} and {A{du{ts)/dt)} are uniformly bounded. We again use 
Lemma 2.5 in [54] to get that Au{ts) Au{t) and A{du{ts)/dt) A{du{t)/dt) as 


Lemma 3.6 Au{t) and A{du{t)/dt) are Bochner integrable on J. 

Proof. We introduce the operator A^ and J™, defined from X into X, given by 

= (^I + ^A^ X 
AF'x = AJ^x for x e X 

Then, for all x and y in X, 

\\J^x-J^y\\ < \\x-y\\ 

\\A^x-A^y\\ < 2ml|x-y|| 

and for x € D{A) 


Now, 


IJA’^zlI < ma:|| (cf. Kato [54]). 


< 


WA'^U^it) - + ||A"‘X”^(t)|| 

2m. — H C 
m 
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Hence is uniformly bounded. Lemma 2.5 of [54] implies that 

Au{t) uniformly on J as m — > oo. Now we can imitate the proof of Lemma 4.6 of 

[54] to get the required result. 

Remark 3 Lemma 3.6 also hold good if the operator A is replaced by the operator 
B. 


We can rewrite equation (3.3.2) as 
dV'^{t) 


dt 


+ AY^{t)+BX^{t-k) = a.e. on J, 


(3.3.6) 


where 


(3.3.7) 


no) = m+aUo + PUi 

f^{t) = f{tf) + auf_i + for t € (t™ i, tf], j = 1,2,- ■■,m. 

Integrating (3.3.6) over (0, t), we get 

V^{t) = t/i - /* AY^(s) ds - f BX'^is - k) ds 
Jo Jo 

+ / r(») ds. 

Jo 

Proof of Theorem 3.1 Prom (3.3.7), for every x* eX*, we have 

< V^{t),x* > = <Ui,x*>- f < AY^{s),x* > ds 

Jo 

- [ <BX”^{s-k),x*>ds 
Jo 

+ Q’r{s),s:’) ds. 

Passing through the limit as m — > oo, using the bounded convergence theorem and 
Lemma 3.5, we get 

<v{t),x*> = <Uux*>-f <Av(s),x*>ds 

Jo 

- f <Bu{s),‘. 

Jo 


,x* > ds 


+ f{s) + au{s) + ^v{s), X * ) ds. 
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The integrand on the right are continuous on J for each fixed x*, hence < v{t),x* > 
is continuously differentiable on J. Using Lemma 3.6 and the boundedness of Au{t), 
Bu{t) and the continuity of /, we deduce that the strong derivative of u(t) exists 
and 

+ Bu(t) = /(t) + au(t) + a.e. on J. (3.3.8) 

Qih 

Since we have 

(jnj 

«(0) = C/o, u(0) = —(0) = Ui and 

at 

, s du, , 

v{t) = -g(t). 

the equation (3.3.8) implies that u{t) is a strong solution to (3.2.2). 

Uniqueness. We may assume without loss of generality (or we can add el to A 
for suitable constant e) that the resolvent set p{—A) contains the origin, i.e., —A 
is invertible. It is easy to show that a strong solution to (3.2.2) is a mild solution 
to (3.2.2). Now we will show the uniqueness of a mild solution to (3.2.2) which will 
imply the uniqueness of a strong solution to (3.2.2). 

Let Ui, i = 1, 2, be two mild solutions to (3.2.2). Let dui/dt — Vi and —Aui = Xi for 
i = 1, 2. Then 

Xiit) = -AUo + iT{t)-I)Ui 

+ [ {T{t - s) - I)G{s, {-Ay^Xiis), Vi{s)) ds, 

Jo 

Vi{t) = T(t)Ui+ f T{t-s)Gis,{-A)~'^Xi{s),Vi{s)) ds, i = l,2. 

Jo 

Hence 

llxi(t) - l2(t)[l + IKW-<'2(i)|| 

< O f [|lll(s)-l2(s)|| + ||!Ji(s) -t;2(s)||] rfs. 

Jo 

Applying Gronwall’s inequality, we get 


xi = X2 and on J. 
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Hence 

u\ = = {-A)~^X2 = M2- 

This completes the proof. The existence and uniqueness result of Theorem 3.1 
guarantees that there exists a unique function u{t, x) satisfying initial and boundary 
conditions (1.5.24) and (1.5.25) such that for a.e. x E (0,1), the function t 
u(t,x) is continuously differentiable in (0,r], and the function t {du/dt){t,x) is 
differentiable a.e. on (0,r) and the equation of motion (1.5.22) is satisfied a.e. on 
(0, T) X (0, 0 for all pairs {Uq, Ui) E W^>p{0, 1) n 1). 


Chapter 4 

Viscoelastic Systems With Finite 
Memory 


4.1 Introduction 

Consider one dimensional motion of a linear non-aging viscoelastic medium (cf. 
Drozdov [56] for more details). The medium occupies the half space a; > 0 initially. 
At time t = 0, the boundary surface x = 0 starts moving in the x-direction. We 
assume that the displacement uo of the boundary is time independent. The body 
forces and surface traction vanish. 


The displacement vector u = vJ'Ei has only one non-zero component v} = u{t,x). 
The tensor e = e^^eiCj of infinitesimal strains also has only one non-zero component 

_ du{t, x) 

^ dx 

Therefore we have 

du{t,x) _ 2 du{t,x) . 

^ “ ax ’ 3 ax ’ ^ 

where e is the first invariant of the tensor e and en is the only non-zero component 

of the deviator e of the strain tensor e. 

The constitutive equation for the non-aging viscoelastic material is given by 


a{t) = 3Keit) 


(4.1.2) 
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s{t) = 2G 


~ [ R{t — s)e{s)ds 
Jo 


+ 2r] 


de{t) 

dt 


(4.1.3) 


where a is the first invariant of the stress tensor < 7 , K and G are volume and shear 
moduli of elasticity, 97 is Newtonian viscosity, R{t) = — Qo(t) is a relaxation kernel, 
Qo{t) is a relaxation measure, here prime denotes the differentiation with respect to 
time. Prom (4.1.1) and (4.1.3), we obtain the following non-zero components of the 
stress tensor a, 

du{s, x) 


,11 


+ 


4 d^u{t,x) 
3^ dtdx ’ 


R{t-s)- 


dx 


ds 


,22 


= 


K 

K 


du{t, x) 
dx 

du{t, x) 
dt 


The equation of motion in the absence of body forces is given by 

d^u 


(4.1.4) 

(4.1.5) 

(4.1.6) 

(4.1.7) 


where p is the material density, V is the Hamilton operator in the initial configura- 
tion and the dot stands for the scalar product. Prom (4.1.6), we obtain the following 
integro-differential equation for u{t,x): 


i’-W~ = + + 

rt 


AG 

3 


[ R{t-s) 
Jo 


d'^u{s, x) 


dx"^ 


ds. 


The initial and boundary conditions are given by 

9u(0, x) 


u(0,x) 
u{t, 0) 


dt 

Uq {t > 0). 


= 0, {x> 0), 


(4.1.8) 

(4.1.9) 
(4.1.10) 


4.2 Abstract Formulation 

With simple modification of the relaxation measure Qo{t) we may write the integro- 
differential equation (4.1.8) as the following abstract integro-differential equation in 
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a Hilbert space H: 

r 

+ -4 aut{t) + g{0)u{t) + / g'{s)u{t — s)ds = 0 (4-2.1) 

Jo J 

where o; > 0 is a constant and A is a linear positive definite, self-adjoint unbounded 
operator on H with the domain D{A) contained in another Hilbert space V, where 
V is densely and continuously embedded in H. On the kernel g we assume the 
following conditions: 

(Gl) g e ( 72 ( 0 , 00 ) n C[0,oo) and g' € L^{0,ooy, 

(G2) g(s) > 0, g'(s) < 0, g" > 0 on (0, oo); 

(G3) the equilibrium elastic modulus g(oo) > 0, which is taken as unity without 
loss of generality; 

(G4) there exist positive constants S, Si and Cg such that 

g"(s)+Sg'(s) > 0, s€(0,oo) 

and 

g"(s) < Cpl 5 '(s)|, sG[si,oo). 

We note that, for instance, the weakly singular kernel of the form 

g-C25 

g'(s) = ci,C 2 > 0, 0 < C 3 < 1, (4.2.2) 

satisfies (G1)-(G4). 

The extensive treatment of linear viscoelastic materials with short as well as long 
memory using variational methods is considered by Glowinski, Lions and Tremolieres 
[31], Duvaut and Lions [32] and others. The particular case of (4.2.1) in which 
A = -A (A is the Laplace operator) and a = 0 has been considered by Dafermos 
[33] subject to the Dirichlet boundary conditions where the history kernel g satisfies 
(G2), (G3) and 

(Gl)' g' e CHIO, oo), s', s" € LH0,oo), 

with s' assumed to be convex and showed that the energy associated with the system, 

E(t) = 


f b'(s)l IKt) - u(i 

Jo 


yds, 
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goes to zero asymptotically. The three parts in E{t) are potential, .kinetic and 
memory energy. 

Results of Dafermos [33] were improved in Dafermos [34] by dropping the convexity 
condition on g'. Day [35] has given some explicit rates of decay of the energy for 
similar problems. Desch and Miller [36] - [37] considered (4.2.1) (with a = 0) under 
the conditions (Gl)\ (G2) and (G3) and zero initial history. The assumption (Gl)^ 
excludes the singular kernels of the type mentioned in (4.2.2). With additional 
smoothness on the kernel and zero initial history, Hannsgen and Wheeler [38], [39] 
established exponential decay of the energy of the system. Fabrizo and Lazzari 
[40] investigated a three dimensional viscoelastic system with memory with the 
conditions (Gl)-(G3) and 

g"{s)+6g'(s) > 0, se(0,oo) 

for some constant 5 > 0. Liu and Zheng [41], [42] investigated (4.2.1) without 
damping, i.e., a = 0, under the conditions (Gl)-(G4). Our aim is to generalize the 
results of Liu and Zheng [41], [42] for strongly damped linear viscoelastic systems 
with memory (i.e., a > 0 in (4.2.1)). 

Let i? be a Hilbert space with the norm || Hij and the inner product (., .)h and let V 
be another Hilbert space with the norm ||| |||v and the inner product < ., . >v such 
that the embedding V H is dense and continuous. Let V* denote the dual space 
of V of all continuous linear functionals on V. Then the triplet {V, H, V*) satisfies 
the continuous embedding 

V^H-^V* (4.2.3) 

where we have identified H* with H itself. Consider a symmetric sesquilinear form 
cr on y satisfying 

|£r(u,u)| < lll^’llk for all u,v€V, (4.2.4) 

a{u,u) > for all ueV, (4.2.5) 

where u and 0^ > 0 are constants. The form a gives rise to a bounded linear 
operator A : V — > V* through the relation 

a{U,V) = y 


for all u, u € V, 


(4.2.6) 
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where < f,v >v,v is the duality paring between f €V* and v £V. The operator 
A : D{A) — > H given by 

D{A) = {u € y : AueH} (4.2.7) 

is called the restriction of A on if and it is a positive definite self-adjoint operator. 
Also, V = D (At) and 

a{u,v) = ^A'^u,A^v^^ for u,veV. (4.2.8) 

We assume that this is the kind of operator we have in (4.2.1). 

The space V can also be equipped with an inner product 

{u,v)v = cr{u,v) for u,v€:V (4.2.9) 

and the inner product (., .)v is equivalent to the inner product < ., . >v. Henceforth, 
on V we shall use the equivalent inner product {.,.)v and the corresponding norm 
||u||y = {u,u)v only. In place of (4.2.4), we shall use 

la-(u,u)| < C'<r|lu||v ||ul|v for all u,v€V, (4.2.10) 


where > 0 is another constant. 

We set W = L‘p{0, oo; V) which is a Hilbert space of all V-valued, square integrable 
functions defined on the measure space ((0,oo), V, l^'|ds) equipped with the norm 

= r\9'{s)\\\w{a)\\lds. (4.2.11) 

Jo 

Let = V X if X W be the product space which is a Hilbert space equipped with 
the norm 


^\\h = ll^llv + 11^111- + llwllwi for 2 : = {u,v,wY eu. (4.2.12) 


We set 


D{A) = ^zen\av + u- g'{s)w{s)dseD{Ay, veV; DsW€W;w{0)=qY 

(4.2.13) 

\ 


Az = 


/ 


V 


-A (oiv + « - /o°° g’{s)w{s)ds) 

V - DsW 


(4.2.14) 
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Then (4.2.1) can be written as an abstract Cauchy problem in U, 

dz 

— ^ Az, t>0, (4.2.15) 

^( 0 ) = zo, 

with the help of the transformations 

V = Ut, w{t,s) = u{t) — u{t — s), t, s€(0, oo) (4.2.16) 

and the initial point taken as 

•^0 = {<!>,, i’, en. (4.2.17) 


4.3 Existence and Uniqueness 

For the case of without damping, i.e., a = 0, it was proved by Liu and Zheng |42] 
that the operator A generates a Co-semigroup of contractions (cf. Pazy [21] for 
definitions). This result was proved with the help of the following lemma stated and 
proved in [42]. 

Lemma 4,1 Suppose that function g satisfies the conditions (Gl)-(G3). Then for 
any w eW with w(0) = 0 and DgW 6 W, we have 

|s5'(s)|-4 0, as s — >• 0, 

l^'(s)| lli£;(s)||y -> 0, as s 0, 

f /(s)||u;(s)i|y ds < 00 , 

Jo 

ly'(s)| llu;(s)l|v->0, as s^oo. 

We prove the following theorem which extends the results of Liu and Zheng [42]. 

Theorem 4.2 Let g satisfy (G1 )-(G4). Then the linear operator A given by (4-2.13) 
and (4-2.14) is the infinitesimal generator of a Co-semigroup of contractions on “H. 

Proof If we show that A is dissipative and the range oi XI — A = TL for some 
A > 0, then from Lumer-Phillips Theorem (cf. Pazy [21] on page 14), it will follow 
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that A is the infinitesimal generator of a Co-semigroup of contractions on H. For 
z 91, we have 


{Az,z)y^ — {v,u)y+^—A^v-\-u — J g'{s)w{s) ,v^ 

-1- (u JD^W^W^y^. 


(4.3.1) 


We note that 


— {Aiav + u 


poo 

/ 5'(s)w(5) 

Jo 


ds^ , = — a + u — J g'{s)w{s) ds, 

= —aa{v, v) — a{u, v) 

poo 

+ / g'{s)cr{w{s),v) ds. 

Jo 

poo 

< —a{u,v)+ I g'{s)a{w{s),v) ds 
Jo 

poo 

= -(«,'yV+ / g'{s){w{s),v)y ds. 

Jo 


Also, since p' < 0 on (0, oo), we have 

poo 

{v,w)yy = / |^'(s)| (u,u;(s)V ds 

Jo 

poo 

= - 9'{s) (v, w{s))y ds. 

Jo 

Hence, from (4.3.1) and (4.3.2), we have 


( 4 . 3 . 2 ) 


{Az, z)-f^ < - (L>sW, 

poo 

= - |/(s)| (■Dau;(s),u;(s))y ds 

Jo 

= I P g\s) ^\\w(s)fy d3 
= r /(B)\\w{s)ry ds 

< 0. (4.3.3) 

Hence the operator A is dissipative on %. Next we show that the range R{A) = %. 
For any F = ^ we want to find a z = (u,u,ty)^ € R such that 

W T jvfhr 
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Az = F which is equivalent to the following system 


—A ^av + u — J g'{s)w{s) ds^ 


V — DcW 


= A, 
= A) 
= A- 


(4.3.4) 

(4.3.5) 

(4.3.6) 


Prom (4.3.4), we have v = fi eV uniquely determined and from (4.3.6), we have 

w{s) = sfi- [ A(t) dr. (4.3.7) 

Jo 

We note that the function A = A from (0? oo) into V is in W since 

poo 

Wfifw = / |5'(s)|llA(s)|P,.ds 

Jo 

= ll/ifv- [f Is'MI *) . 

Therefore from (4.3.6) we have that DgW = fi + fz € W. Also, we have u;(0) = 0. 
Now as in [42], we show that w €.W. For any T, e > 0, we use (G4) and Cauchy- 
Schwartz inequality to get 

J l5'(s)l lh(s)llv 

= ]9\T)\\w{T)\\l--^9\e)\\w{e)\\^ 

2 /*^ 

~ S J 

+ ^r\g'(^)\\\DM\lds ( 4 . 3 . 8 ) 

Prom (4.3.8), we have 

J i5^(s)i im^)iiv 

+ -fimWOMlds. ( 4 . 3 . 9 ) 
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Taking the limits as € -> 0 and T ->• oo in (4.3.9) using Lemma 4.1, we obtain 

o 4 

Ikliw < |p'(5)| ||Dsiu|l^ w ds 

= ^ \\Dsw\\l, (4.3.10) 

which implies that w e W. Now, since A : D{A) — > H is onto, let h € D{A) be 
such that Ah = / 2 . Then 


u 


h — av + 


f 


g'is) 


sfl 


f Mr) 

Jo 


dr 


ds 


(4.3.11) 


solves (4.3.5) uniquely. Thus, for any given F there is a unique z 6 D{A) such 
that Az = F. Now, we show that there exists a positive constant Co such that 


\ z\\h < C-oIIFIIk. 


(4.3.12) 


From (4.3.5), we have 

- (^A ^Q!U + u - y* g'{s)w{s) ds'^ , = (/a, u)h 

which is equivalent to 

-cr(QV + u- J g'{s)w{s) ds, uj = (/a, u)h 
and after rearranging the terms leads to 

poo 

—{f 2 ,u)H — ota{v,u)+ / g'{s)a{w{s),u) ds = cr(u,u). 

Jo 

Prom (4.2.4), (4.2.5) and (4.3.13), we obtain 

ll^^llv = <r{u,u) 

< 11/2|UI|w||h + 0!C'(rllu||v||«|ly 

/ oo 

l5'(s)l l|wj(s)l|vl|w||y ds 

Cff f l^'(s)lh{s)llvll^itlv ds 
Jo 


(4.3.13) 


(4.3.14) 




(4.3.15) 
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where Cvh is a positive constant for the continuous embedding V ^ H such that 

lltillff < C'vffllujly. 

We use the inequality 2a6 Ka^ + b"^ to get 



Similarly, we have 


aC'<j(|u||y|lu||y < 




and 


op2 

a||u;(s)||yllu||y < % 


^Nlv+^IN 




(4.3.16) 


(4.3.17) 


+ 




-1 


(4.3.18) 


The inequalities (4.3.18), (4.3.10) and the equation (4.3.6) together with the fact 
that V = f I imply that 


roo 

C. / |5'(S)1 i|t^;(s)llvllu|ly ds 

Jo 


ZCl 
“ 2 

< ^ 
- (52 

12CI 


|/(s)|ds| + 


sMlw + gll^^llv 


< 

< 


- (52 

12CI 

(52 


+ gll^llv- 


([y 

(^ |5'(«)I <(s) ll-D 

(j" |/(s)| ds'j [llAlll, + IIAII^] + l||«i|“v 
(jf“ U'(s)| ds) [(j(“ |s'(s)| ds'j ll/ill?. + II/, II 


(4.3.19) 


Combining the inequalities (4.3.15), (4.3.16) and (4.3.19), we obtain 

2l 




Za^Cl + 


24(72 

^2 




liAII 




SCjallAlpH + ^ (jf 
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Also, we have 

li^llH = ll/i|llr<<^VHll/illv (4.3.21) 

and 

llwllii. < |llAw|||^ 

~ ^ Uo ■ (4.3.22) 

Prom the inequalities (4.3.20)-(4.3.22), we have 

IWI« < C$H + 3a=Cj+|(^°°|j'(s)|&)+^(^”|ff'(s)|<is)' ll/illf, 

+ ^ (l + 3C^) (^ lff^(5)l tis) ll/all^. (4.3.23) 

The inequality (4.3.23) implies that there exists a positive constant Co independent 
of z and F such that 

ll^lk < CollFll^. (4.3.24) 

Thus A is invertible and A~^ is in the space B{'H) of bounded linear operator from 
H into itself. Now, from Theorem 5.2-A in Taylor {57] on page 260, it follows that 
for any A > ||>1“^|1 b(h)) the operator XI - A~^ is invertible and (XI — is in 

B{7{). Therefore the operator XI - A is invertible and belongs to B('H) since 

XI-A = -A{I-XA-^). (4.3.25) 

The Lumer-Phillips Theorem now gives the required result. This completes the 
proof of the theorem. 

The following result due to Hille [58] (cf. also, Pazy [21] on page 102) for an 
abstract Cauchy problem of the type (4.2.15) ensures the existence of a unique 
continuously differentiable solution. 

Theorem 4.3 Let A be a densely defined linear operator on U with nonempty re- 
solvent set p(.4.). The abstract Cauchy problem (4-2. 15) has a unique solution in 
C^([0, oo); “H) for every Zq G 11(A) if and only if A is the infinitesimal generator of 
a Co-semigroup onM. 
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4.4 The Exponential Stability 

In the present section we prove that the semigroup S{t) of contractions on % gen- 
erated by A. is exponentially stable. This result extends the corresponding result of 
Liu and Zheng [42]. 

To establish the exponential stability of the semigroup S{t), we require the following 
theorem (cf. Theorem 1.3.2 in [42] on page 4 ). 

Theorem 4.4 Let S{t), t > 0 be the semigroup of contractions on % generated by 
A. Then S{t) is exponentially stable if and only if 

p{A) 2 {iPJeR} = iR (4.4.1) 

and 

limsup|^l_,«,||{W-.4)"Nls(«) < oo. (4.4.2) 

where R is the set of reals. 

We have the following result. 

Theorem 4.5 The semigroup S{t), t > 0 of contractions on % generated by A is 
exponentially stable, i.e., there exist positive constants M and 7 such that 

|| 5 (t)|| < for all t>0. 

Proof As we have already shown that A is invertible and A~^ belongs to 

the operator ipA~^ — / is invertible for G R with \P\ < Hence the 

operator 

01 -A = A{0A~^ - 1) 

is invertible for j^j < 2 nd \\{0I — .A)~^||b(k) is a continuous function of 

P in the interval (-M“^ 11 b(-«)> ll-^“^lli(«))- 
If 

sup{||(i/3/ - .4 )"^||b(«) 1 1^1 < = M < 00 , 

then by contraction mapping theorem, the operator 

01- A = {i^I-A)iI + i0-0o)i0oI-A)-'^) 
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is invertible for \Po\ < \\A and - /?o| < (1/M). For \0q\ very close to 

{p I 1^1 < + (l/M)} c p{A) 

and \\{iPI - >A)“^1|b(k) is a continuous function of p in the interval 

- (1/M), + (1/M)). 

Therefore if (4.4.1) is not satisfied, then there exists cj € R with < |a^l < 

00 such that {0 | \p\ < |u;|} C p{A) and sup{ll(ij0/-^)“^l|B(-«) | |;0l < \uj\} = oo. It 
follows that there exists a sequence € R with ->■ w, j/?„| < |a;| and a sequence 
of complex vector functions = {un,Vn,Wn) G D{A) such that 

ll^nlll, = llt^nlPv + lbfff+INlPw = 1 


for all n and 


(iPnl - A)zn -^0 in "H as n oo. 


(4.4.3) 

The convergence in (4.4.3) is equivalent to 



iPn^n ^ i^ 

V, 

(4.4.4) 

ipnVn + A(avn + Un- f /(s)u;„(s) ds J 0, in 

H, 

(4.4.5) 

\ < 1 / 0 ' 

iPnWn -Vn + Ds^n ->• Oj 

W. 

(4.4.6) 

As in (4.3.1), taking the inner product of Zn with (4.4.3) in H and collecting the 

real parts, we obtain 



Re {Azn,Zn)-u = /(5)|K(s)|Ivds^0. 

(4.4.7) 


Prom (4.4.7) and (G4), it follows that t;„ ->■ 0 in JT and u;„ -> 0 in W. Hence 

IKiPv 1- 

Now, taiiag the inner product of % with (4.4.4) in H, and using the fact that 
u„ -> 0 in H, we get 


iPn{Un,'f}n)H 


(4.4.9) 



Chapter 5 

Regularity of Solutions 


5.1 Introduction 


Consider the following initial value problem for the semi linear hyperbolic integrod- 
ifferential equation in the n-dimensional Euclidean space R", 

+ k{t-s)g{s,u{x,s), ^^^^^'^ )ds, t > fo, a: e R", (5.1.1) 


u{x, to) 


Ui{x), — - 1 " - ° -- = U 2 {x), X € R”, 


where A is the n-dimensional Laplacian, / and g are smooth nonlinear functions 
and fc is a locally p-integrable function for 1 < p < oo. Problem (5.1.1) is equivalent 
to the first order system 


d f u{x,t) \ 
^ \ v{x,t) j 



uix,t) \ / 0 \ 

) \f{t,uix,t),v{x,t)) J 

^ ] ds, t > to, X G BP 

g{$,u{x,s),v{x,s)) J 


and 


/ u{x,to) ^ 

\ v(x, to) J V U 2 ix) 


x€R". 
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Thus, for some appropriate choice of a Banach space X of functions, problem (5.1.1) 
can be treated as a particular case of the following integrodifferential equation in 


du{t) 


+Au(t) = f{t,u{t)) + k{t- s)g{s,u{s)) ds, £ > to, (5.1.2) 
u(£o) = uq. 


where A, defined from the domain D{A) C X into X, is a linear operator, the 
nonlinear maps / and g are defined from Jj- x X into X, with Jt being the closure 
of the interval [tojT), £o < T < oo, and the real- valued map k is defined on Jt- 
For (5.1.1), the right choice of the space is the Hilbert space H = JJ^(R") x L^(R”). 


The operator —A is associated with the matrix of operators 


0 I 
A 0 


and it is 


known that this matrix of operators is the infinitesimal generator of a Co semigroup 
of bounded linear operators in H (cf. Pazy [21]). 


Initial studies on problem (5.1.2) for the case in which fc = 0 was considered by 
Segal [59]. Different forms of solutions in this particular case have been considered 
by Martin [18] and Pazy [21] . An example in which the mild solution is not a strong 
solution is given by Webb [60]. Further modifications of results of Segal [59] have 
been considered by Kato [54], Martin [18] and Pazy [21]. 

Our aim here is to establish the existence and uniqueness of regular solutions 
to (5.1.2) which will in turn guarantee the existence and uniqueness of regular solu- 
tions to (5.1.1). We first prove the local existence and uniqueness of a mild solution 
to (5.1.2) under the assumptions that -A generates a Co semigroup of bounded 
linear operators in X, fit,u) and g{t,u) are continuous in £ and satisfy certain 
local Lipschitz conditions in u with Lipschitz constants depending on £ and ||u|]x 
and k e LP{0,T), 1 < p < oo. Under the assumption of continuity on k on [0,T) 
we show that if the local mild solution exists on Jti = [£o,^i], to < ti < T, then 
it can be extended further. Also, we show that either the solution exists on the 
whole Jt or there is a maximal interval of existence [£o, tmax)^ ^0 ^rrmx ^ 3»Ild 
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ljti(4)||A- = oo. Furthermore, we show that (5.1.2) has a classical (continu- 
ously diflferentiable) solution provided / and g are continuously differentiable from 
[to,T) X X into X. 

For the various regularity results in the case when the operator —A in (5.1.2) gener- 
ates an analytic semigroup, we refer to Rankin [61] and Bahuguna [49] and papers 
cited therein. Note that for the applications to hyperbolic problems, we cannot 
assume that —A generates an analytic semigroup. 

The results presented here can be generalized for the case of time dependent oper- 
ator A{t) using the theory of evolution systems. For simplicity of presentation, we 
restrict ourselves to the case A{t) = A. 


5.2 Main results 


Let X be a Banach space and let ||.||x denote the norm in X. Let to > 0 and let 
to<T <oo. Let denote the closure of the interval [to, v) for any to <v<co. Let 
C{J,j : X) be the space of all continuous functions from Jr, into X. For a compact 
interval J„, let C{Jr , : X) denote the Banach space of all continuous functions from 
Jr, into X endowed with the supremum norm 


1|/||C(J,:X) = sup ||/(t)||jf, f e C{Jr, : X). 

t^Jr\ 


It is known that if —A is the infinitesimal generator of a Cq semigroup T(t), t > 0 
and (5.1.2) has a classical solution or a strong solution u on Jt then u satisfies the 

integral equation 


u{t) = T{t - to)uo + j^T(t- s)[fis, u{s)) 

-t- f k{s- r)g{r,u(r)) dr] ds, t G Jt- (5.2.1) 

Jto 

Clearly, u given by (5.2.1) is in C{Jt : X). Hence vre define the mUd solution 
to (5.1.2) as follows. 
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Definition 5.1 A function u € C{Jt : X) satisfying the integral equation (5.2.1) 
is called the mild solution to (5.1.2) on Jr. 

We have the following result for the mild solutions to (5.1.2). 

Theorem 5.1 Let —A be the infinitesimal generator of a Cq semigroup {T{t),t > 

0}, in X. Let f ,g : Jt "X- X -i- X be continuous in t on Jx and satisfy the following 
conditions. 

(HI) For every compact subinterval Jx^ of Jx and allx,y e Bii(0,A') = {z e X : 
< R} there exists a constant Li{To,R) such that 

\\f{t,^)-f{t,y)\\x < L,(To,R)\\x-y\\x. 

(H2) For almost every t E Jx o,nd all x,y E Bji{0,X) there exists a nonnegative 
function R) E L^^^{Jx), 1 < p < oo such that 

|lp(t,x) -p(t,?/)llx < Li(t,i?)l|x-y|lx. 

Further, suppose that 

(HS) The real-valued map k is in L^(0,r) where 1 < q < oo is such that 1/p+l/q = 

1 . 

Then there exists a unique mild solution u to (5.1.2) on a compact subinterval Jt^ 
of Jx- If k is continuous on [0,T), then either the mild solution u can be extended 
to a unique mild solution to (5.1.2) on Jx or else there exists at^ax, to < tmax < T, 
such that (5.1.2) has a unique mild solution u on [to,tmax) ond 

Urn |lu(f)|| = oo. 

tjtmax 

Under stronger assumptions of differentiability on the nonlinear functions / and g, 
we have the following regularity result establishing the existence and uniqueness of 
classical solution to (5.1.2). 

Theorem 5.2 Let A be the infinitesimal generator of a Cq semigroup {T{t),t > 0} 
on X. If f,g : Jx 'x X X are continuously differentiable from Jx x X into X 
and k is continuous on Jx, then the mild solution u to (5.1.2) with uq E D{A) is a 
unique classical solution to (5:1.2) on Jx- 
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5.3 Proofs of Main Results 


Proof of Theorem 5.1 : We fix an arbitrary f,to<t <T. We first prove the 
existence and uniqueness of a mild solution to (5.1.2) under the assumptions that 
Li{f,R) and L 2 {t,R) are independent of i?, i.e., Li{f,R) = Li{f) and L 2 {t,R) = 
L 2 {t) for all R> 0. Let 

M{t) = sup ||r(t)||B(x), 

Mo(t) = 1 + M{t){Li(f) + ll^l|LP(Jr)ll-^2||L«(JT)) 

Now, define a map F : C{Jf : X) C{Jf : X) by 

{Fu){t) = T(t-to)uo+ [*T{t-s)[f{sMs)) 

Jto 

+ / k{s - T)g{T,u{T)) dr] ds. ( 5 . 3 . 1 ) 

Jto 

Now, for u,v e C{Jf: X), we have 

||(Fu)(t) - (Fu)(t)|U < Mo(T)(t - to)\\u - 


Repeated applications of the above inequality gives us 

\\{F^u){t) - (F"u)(t)|U < ^[Mo(f)(i - to)]"||u - v\\c{jf,x). 
Therefore we have 


1 


\\F^u-F^v\\cijf:X) < -:[Mo{T)'IT\\u-v\\c{j^..x). 


m 


(5.3.2) 


(5.3.3) 


(6.3.4) 


Prom the Banach contraction principle, it follows that there is a unique fixed point 
u € C{Jf : X) of the map F. This u satisfies the integral equation (5.2.1) and 
hence it is a unique mild solution to (5.1.2) on Jf. Moreover,’ if v is the mild 
solution to (5.1.2) with u(to) = on then 

||u(t) — u(t)|lx’ < Mo(T')||uo - + M)(T) f llu - v^c{^j,-.x) ds. (5.3.5) 

J to 

Thus, for every r) & Jf, 'we have 

11^(77) — u(77)}ix ^ MoiT)\\u(^ — uolU + M)(T) f It^ - ds (5.3.6) 

J to 
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Taking the supremum over Jj, we get 

||u - < ^o(T)|l«o-i;o|U + Mo(r) /* ||u - uJ|c(j,:X) (5.3.7) 

JtQ 

Applying Gronwall’s inequality and taking the supremum over Jf, we obtain 

- 'y||c(jj.:X) < ■Wo(r)exp|Mo(f)f}||uo -uo)jU. (5.3.8) 

The inequality (5.3.8) proves the uniqueness and continuous dependence of a mild 
solution to (5.1.2) on the initial data on Jf. 

Now, under the general assumptions (HI) and (H2) on / and g we show that there 
exist a unique mild solution u to (5.1.2) on Jt^ for some ti, to < h < to + 1. 
Furthermore, under the assumption of continuity on k on Jt we show that either the 
solution can be extended on to the interval Jt or else the solution can be extended 
on to the maximal interval of existence [to,^max), to < tmax < T, and 


Let 


lim I[u(t)||A: = oo. 

Mima® 


R{to) — 2M(to + l)||uo|U) 

Ni{to) = sup l|/(s,0)||x, 

to<t<to+l 

N 2 {to) = sup ||p(s,0)|U, (5.3.9) 

C{to) = R{to)Liito + 1) R{to)) + Ni{to) 

+ ■R(io)l!fc||LP(Jf)ll-^2(-5J?(io)!U«(Jii) +-^2(io)Plliy'{Ji^)- 

We set 

5(k,uo) = (5.3.10)' 

and let ti = to + ^{to, uq)- 

Now, we show that the mapping F given by (5.3.1) maps : X)) into 

itself. For u € BR(to){0,C{Jt, : X)), we have 

l|F(u)(t)l|x < M(to + l)|MU+i^(*o + l) f [ll/(s. ■«(«)) -/(s>0)lk + ll/(s,0)!l 

V to 
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+ f l^(s -t)\[ ll5(T,M(r)) -5(r,0)j|x + 1|5 (t,0)|U] dr] ds 

< M(to + 1)[ l|tio|U + C'(to)](i - to) 

< 2M(to + l)I|uo|U = i2(£o)- (5.3.11) 

Also, for u,v €. C(Jtj : X)), we have 

l|Fu - < M(£o + l)C(to)(£i - to)\\u - y||c( 7 *,:X)- (5.3.12) 

As before, the above inequality after repeated applications gives a unique fixed point 
u in Bij(to)(0, C( Jti : X)) of F as before. Now, let 


fh 

r{t) = / k{t-T)g{T,u{T)) dr, 

Jto 

= f{t,u)+r(t). 


(5.3.13) 

(5.3.14) 


Consider the initial value problem 
dw{t) 


dt 


+ Aw{t) = fi{t,w{t)) + f k{t — s)g{s,w{$)) ds, i> £i, (5.3.15) 

Jti 

w{ti) = u{ti), 


Since fi is also continuous and satisfies (HI) with the same Lipschitz constant, 
therefore if ti < T then proceeding as before, we can show that there exists a unique 
mild solution w to (5.3.15) on [ti,t 2 ], £2 = + d{ti, j|u(£ijlA-) where 5{ti, |ju{£i)|lx) 

is defined as 5{to, |luo|lx) by replacing £0 and l|wolU hy £1 and ||u(£i)Hx, respectively. 
Then the function ui defined on Jt^ in to X given by 


^'^^~\ «;(£), £€[£ 1 , £ 2 ] ’ 


is the mild solution to (5.1.2) on Jtj. We need to check this only for [£i,£ 2 ]- If 
£ G [£i, £ 2 ] then 

u(£) = w{t) = T{t - ti)u{ti) + f Tit- s)[/i(s, u;(s)) 

Jtl 


+ 


f kis- r)g{'i 
Jtl 


,w{t)) dr] ds 
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+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


T{t ii)[T(ti — to)tio + f T{ti — s)[f{s,u{s)) 

/ k{s - T)g(T, u(r)) dr] ds] 

^ to 

f T{t- s)[fi{s,w{s))+ f k{s - T)g{T,‘w{T)) dr] ds 

Jti Jti 

rh 

T{t-to)uo+ - s)[/(s,ui(s)) 

Jto 

/ k{s - T)g{T, Ui(r)) dr] ds 

Jto 

[‘nt-s)if{sM‘)) 

Jti 

r(s) + k{s- r)^(r, tii(r)) dr] ds 

Jti 

T{t-tQ)uo+ [ T{t — s)f{s^u{s))d$ 

Jto 

rh rs 

I I T(t — s)k{s — r)g(T, Ui(r)) dr ds 
Jto J to 

ft fti 

I / T{t — s)k{s — r)g{T^ Ui(t) dr ds 

Jti Jto 

n T{t - s)k{s - r)g{r, Ui{t)) dr ds 

X 

T{t-tQ)uo + [ T{t- s)[f{s,uiis)) 

Jto 

/ k{s - T)g{T, Ui{t)) dr] ds. 

Jto 


Thus, either there exists a unique mild solution u to (5.1.2) on the whole of Jt or else 
there exists a unique mild solution on the maximal interval of existence \tQ,tmax)i 
to < tmax < T. Now, if limttt„„. lk(i)lU < oo, then as before we can extend the mild 
solution beyond tmax which would contradict the definition of tmax- This completes 
the proof of the theorem. 

The proof of Theorem 5.1 can be modified to get the following result. 


Corollary 5.3 Let A, f, g and k be as in Theorem 5.1. Let r € C{Jf : X), Then 
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the integral equation 

w{t) = r{t) + f T{t — s){f{s,w(s)) + f k{s — T)g{T,w{T)) dr] ds, t^Jf 
"'*0 Jto 

has a unique solution in C{Jf ; X). 

Proof of Theorem 5.2 : If / and g are continuously differentiable from Jy x X 
into X then for any compact subinterval Jf oi J^, f and g are continuous in t on 
Jf and satisfy (Hi) and (H2). Therefore (5.1.2) has a unique mild solution u on Jf. 
We shall show that u is continuously differentiable on Jf. Let 


B,{t) = -fit,u), 
Blit) = -^gityU), 


(5.3.16) 

(5.3.17) 


(5.3.18) 


r{t) = -AT{t-to)uo + T{t-to)f{to,Uo) 

+ T{t- s)k{s - to)g{to, uo) d$ 

Jto 

d 

+ k(s- r)—g{T, u{t)) dr] ds. 

Jto 

Consider the integral equation 

w{t) - r{t)+ f T{t - s)[Bi(s)w{s) 

Jto 

+ k{s- t)B 2 {t)w{t) dr] ds. 

Jto 

The assumptions on / and g imply r is continuous on Jf and 5i(t, u) = Bi{t)u are 
continuous in t from Jf into X and uniformly Lipschitz continuous in u. Prom 
Corollary 5.3, it follows that (5.3.19) has a unique mild solution w on Jf. Now, 
from the assumptions on / and g we have 

/(s,u(s + /i))-/(s,u(s)) = Bi{s){u{s + h) - uis)) 

+ wi(s, h). 


(5.3.19) 


(5.3.20) 
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g{T,u{T + h)) - g{T, u{r)) = B 2 (r)(w(r + h) - u(r)) 

+ oj2{T,h), (5.3.21) 

f{s + h,u{s + h))- f{s,u{s + h)) = —f{s,u{s + h)).h 

+ ujz{s,h), (5.3.22) 

Q 

g{T + h,u{T + h)) - g{T,u{T + h)) = —g{T,u{T + h)).h 

+ ujA{s,h), (5.3.23) 

where /f^llt!;i(s, h) j|x 0 as /i -> 0 uniformly on Jf for i = 1, 2, 3, 4. Let 

Wh{t) = + (5.3.24) 

n 


Then 


+ h — to)uo — T(t — to)uo) + AT{t — to)uo 

+ T{t + h-$)[f{s,u{s)) + k{s-T)g{T,u{'^))dT]ds 

- T{t - to)f{to, uo) - f T{t - s)k{s - to)g{to, uq) ds 

Jto -1 

i [ T{t + h- s){f{s, u(s)) + [ k{s- T)g{T, u{t)) dr] ds 
^ Uto+h JUi 

- r T{t - s)[f{s, u{s)) + [ k{s- T)g{T, u{r)) dr] ds 

Jto Jto 

- jyit-^)[§;fMs)) + l k{s-T)^g{T,u{r))dT]ds 

_ / T{t - s)[Bi{s)w{s) + kis-T)B2{r)w{T) dr]ds. ( 5 . 3 . 25 ) 

Jto Jto 

Now putting s = 7] + h and then replacing rj by s, we get 

T(t + h- s)[f{s, u(s)) + f k{s- T)g{r, u{t)) dr] ds 
Jto+h J^ 

= f Tit - r))[f{v + h, uin + h)) 

Jto 

prj+h 

+ / kir) + h-T)g{T,uiT))dT]drt 

Jto 
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— f T{t- s)[f{s + h,u{s + h)) 

•J to 
fs-hh 

+ k{s + h- r) 5 (r, w(r)) dr] ds 
Jta 


(5.3.26) 


Again, in the inner integral on the right of (5.3.26) putting r = 7 + h and then 
replacing 7 by r, we get 

ps-hh 

I T{t — s)[/(s + h, u{s + h))+ k{s + h — T)g{T, u{t)) dr] ds 

J to V to 

= I T{t - s)[f{s + h,u{s + h)) + f k{s - + h,u{'y + h)) d'y] ds 

J to Jto^h 

= f T{t — s)[f{s + h,u{s + h)) + f k(s — T)g{T + h,u{T + h)) dr] ds 

J to Jto—h 

= j T{t — s)[f{s + h,u{s + h)) + f k{s — T)g{T + h,u{T + h)) dr] ds 

J to J to 

— f f T{t — s)k{s — T)g{T + h,u{T + h)) dr ds (5.3.27) 

Jto Jto 


Now, using (5.3.27) and (5.3.20)-(5.3.23) in (5.3.25) and readjusting the terms, we 
get 


Wh{t) 


to)uo 


= {T{t + h — to)uo — T{t — to)tio) + AT{t — 
r 1 fto+h 

+ T{t + h- s)[f{s, u{s)) 

+ f k{s - r)g{r, u{t)) dr] ds - T{t - io)/(io, ^o) 

Jto J 

^ f T{t-s) 0Ji{s,h) + uj2{s,h)+ I k{s-r){w3{r,h)+Wi{r,h)} dr ds 

hJto L Jto J 

- s) [|^/(s> + h)- ^/(s, u(s))| 

+ ^ k(s-r)S^^gir,u{r + h))--^g{r,u{r))^ dr ds 

_ jT _ 5) jf g(r + h, u{r + h)) dr + k{s - to)g{to, uo) 

+ r T{t - s)[Biis)wr{s) + [ k{s- r)B2iT)wH{r) dr] ds. 

Jto Jto 


(5.3.28) 
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Since the norms in of all but the term in the last line tend to zero bs h —¥ 0, we 
have 

< e{h)+C{f) [ |lu;/,|ic(j,:x) (5.3.29) 

Jto 

where €{h) 0 as h — >• 0 and 

C(f) = max{||r(i - s)||s(;t)[IIBi(s)ilsW + Pllwwl|B2(s)l|sm] : s e Jf}. 
Applying Gronwall’s inequality in (5.3.29), we obtain 

\H\\ciJr.x) < €{h)exp{Cif)t}. (5.3.30) 


Therefore ||u;ft(t)j|x ^ 0 as h — >■ 0. Hence u is differentiable on Jf and its derivative 
is w on Jf. Since w E C(Jf : X), u E C^(Jf : X). Now, assumptions on / and 
g and u E C^{Jf : X) imply that the maps t ^ f{s.u{s)) and 1 1-4 g{s,u{s)) are 
continuously differentiable on Jf, it follows that 


v{t) = T{t-to)uo+ f T{t-s)[f{s,u{s))-h f k{s -T)g{T,u{T)) dT]ds 

Jto Jto 


is a unique classical solution to 

^ + Av{t) = fit, uit)) + f k{t - s)gis, uis)) ds. (5.3.31) 

dt Jto 

on Jf. By definition, u is a mild solution to (5.3.31) on Jf. By uniqueness of mild 
solutions to (5.3.31), we have u = v on Jf. Thus u satisfies (5.3.31) and therefore 
u is a unique classical solution to (5.1.2) on Jf. Since to <T <T, arbitrary, u is a 
unique classical solution to (5.1.2) on Jy. This completes the proof. 



Chapter 6 


Approximation of Solutions 


6.1 Introduction 

In this chapter we are interested in the Faedo-Galerkin approximation of solutions 
to the following integro-differential equation considered in a separable Hilbert space 
11 , (-,•)), 

u'{t) + Au{t) = /(t,w(t))+/' a{t - s)g{s,u{s))ds, t>0, (6.1.1) 

Jo 

u(0) = (j). 

In (6.1.1), the linear operator A satisfies the following hypothesis. 

(Hi) A is a closed, positive definite, self-adjoint linear operator from the domain 
D{A) C H into H such that D{A) is dense in H, A has the pure point spectrum 


0 < Ao < Ai < As < • • • 


and a corresponding complete orthonormal system of eigenfunction {uj}, i.e., 

Aui = XiUi and {ui,Uj) = 6ij, 
where % = 1 if i— j and zero otherwise. 

If (Hi) is satisfied then —A generates an analytic semigroup in H which we denote 
by e~*^, t >0. It follows that the fractional powers A“ of A for 0 < a < 1 are well 
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defined from D(A“) C H into H (cf. Pazy [21], page 69-75). D{A-) is a Banach 
space endowed with the norm 

IMU = ||A^a:l|, xeDiA^). 

On the kernel a we assume the following condition. 

(H 2 ) a € oo) for some 1 < p < oo. The nonlinear functions / and g are 

assumed to satisfy the following hypotheses. 

(H3) The map / is defined from [0, oo) x Z)(A“) into H and there exists a nonde- 
creasing function Fr from [0,oo) into [ 0 ,oo) depending on > 0 such that 

\\f{t,u)\\ < FR{t), 

\\f{t,u^)-f{t,U2)\\ < F^(f)llni-«2|U 

for all {t,u) and in [0,oo) X Br{D{A°‘),(I>), where BR{Z,Z(i) = {z € 

Z \\\z — Zo\\z < R} for any Banach space Z with the norm || jj^. 

(H4) The map g is defined from [0,oo) x D{A°‘) into H and there exists a non- 
negative function Gr G Lf^g( 0 ,oo) depending on jR > 0 , where 1 < 5 < oo, 
( 1 /p) -1- {1/q) = 1 , such that 

||p(i,M)l| < GR{t), 

l|p(i>Wi) -p(i,U 2 )ll < GR{t)\\Ui-U2\\a- 

for a.e. t G [0,oo) and all u, ui and U 2 in BR{D{A°‘),(j>). Initial studies concerning 
existence, uniqueness and finite-time blow-up of solutions to ( 6 . 1 . 1 ) for the following 
special case of ( 6 . 1 . 1 ), 

u'{t) + Au{t) = h{u{t)), t>0, ( 6 . 1 . 2 ) 

u(0) = <l>, 

have been considered by Segal [59], Murakami [67] and Heinz and von Wahl [ 68 ]. 
Bazley [69, 70] has considered the following semi linear wave equation 

u"{t) + Au(t) = h{u{t)), t>0, 
w( 0 ) = 0 , 

( 0 ) = 


(6.1.3) 
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and has established the uniform convergence of approximations of solutions to (6.1.3) 
using the existence results of Heinz and von Wahl [68]. Goethel [71] has proved the 
convergence of approximations of solutions to (6.1.2) but assumed h to be defined 
on the whole of H. Based on the ideas of Bazley [69], [70], Miletta [72] has proved 
the convergence of approximations to solutions of (6.1.2). The existence, uniqueness 
and continuation of classical solutions to (6.1.1) is considered by Bahuguna [49]. In 
the present work, we use the ideas of Miletta [72] and Bahuguna [49] to establish 
the convergence of Faedo-Galerkin approximations of the solutions to (6.1.1) under 
more general conditions on the nonlinear maps / and g. 

The plan of the chapter is as follows. In the second section, we consider approximate 
integral equations and establish the existence and uniqueness of solutions to these 
approximate integral equations using Banach contraction principle. In the third 
section, we prove the convergence of the approximate solutions to a solution to the 
integral equation associated with (6.1.1). We further show, in this section, that the 
solution can be extended to the maximal interval of existence and it is unique. In 
the fourth section, we consider the Faedo-Galerkin approximations and prove some 
results concerning the convergence of such approximations. Finally, in the fifth 
section, we cite some applications of the results obtained in the earlier sections. 

6.2 Approximate Equations and Solutions 

We shall continue to use the notations introduced in the earlier section. The exis- 
tence of solutions to (6.1.1) is closely associated with the following integral equation 

u{t) = e-^^(f)+ [ e-(*“")^[/(s,u(s)) 

Jo 

+ f a(s - T)g{T,u{r)) dr] ds, t > 0. (6.2.1) 

Jo 

In this section we shall consider an approximate integral equation to (6.2.1) and 
establish the existence and uniqueness of solutions to the approximate integral 
equation. By a solution u to (6.2.1) on [0,r], 0 < T < oo, we mean a function 
u e X.(T) for some 0 < a < 1 satisfying (6.2.1) where X^iT) is the Banach space 
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C7([0, T], £)(i4®)) of all continuous functions form [0,r] into endowed with 

the supremum norm 

ll“IU„(r) = sup ||A“u(t)||. 

0<t<T 

By a solution u to (6.2.1) on [0,r), 0 < T < oo, we mean a function u such that 
u e Xa{T) for some 0 < a < 1 satisfying (6.2.1) on [0,T] for every 0 < T < T. 


Since -A generates the analytic semigroup t > 0, we may add cl to -A for 
some constant c, if necessary, and assume without loss of generality that ||e''*^|| < M 
and that —A is invertible. Furthermore, it follows that commutes with and 
there exists a constant > 0 depending on a such that 


t>0. (6.2.2) 

Let 0 < To < oo be arbitrarily fixed and 

L{R) = (1 + i?)(FR(To) + ||a||j:,p(o,To)||<Jii||j[,«(o,To))- (6.2.3) 

Let 0 < T < To to be such that 

T < min |to, 

and 

sup (6.2.5) 

o<t<r ^ 

Let Hn denote the finite dimensional subspace of Hilbert space H spanned by 
{uo,Uu and let : H — y Hn he the corresponding projection operator 

for n = 0, 1, 2, • • •. We define 


|(1 - a)(i(fi)CJ 


-1 


1— a 


(6.2.4) 


fn, 9n : [0, T] X Xa{T) y H, 

fn{t,u) = 9n{t,u) = g{t,PMt)). 

We set ^(t) = ^ for t € [0, T] and define a map Sn on Bjj(Xci.(T),^) as follows. 

{SnU){t) = [ e-(*-*>^[/n(s,w) 

Jo 

+ f a{s- T)gnir, u) dr] ds. (6.2.6) 
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Proposition 6.1 Let (Hi)— (H4) hold. Then there exists a unique € BR{Xa{T), 
such that SnUfi = Un for each n = 0, 1, 2, • • •, i.e., satisfies the approximate inte- 
gral equation 

Mt) = f 

Jo 

+ / a(s - r)p„(r, u) dr] ds. (6.2.7) 

Jo 

Proof: We claim that Sn : Br{Xc,{T),$) BR{Xa{T),^). For this, we need 

to show first that the map 1 1 -> (5„«)(t) is continuous from [0,r] into Z)(A“) with 
respect to || |la norm. For ti,t 2 e fO,T] with ti < h, we have 


||(5„M)(t2) - iSnu){h)\U < l|(e-‘^^ - 


+ 

r ||e-<*-)'‘.4“|l 1 ||/„(s,u)|l 



Jtl 

r 


+ 

/ \a{s-T)\\\gn{T,u)\\dr]ds 

V 0 



rti 


+ 

Jo 

^^)ll 

+ 

/ |a(s-r)| Il5n('r,«)ll dr] ds. 

Jo 

(6.2.8) 


The first integral can be estimated as 





+ [ ||s„(T,ti)l| tir] ds 

Jo 


< (6.2.9) 

To estimate the second integral, we use part (d) of Theorem 6.13 in Pazy [21] on 
page 74 which states that for 0 < /0 < 1 and x G D{A^), 


11(6-“ - /)ill < Ceti>\\x\\f. 

If 0 < ^ < 1 is such that 0 < a + ^ < 1, then A'^y e D(A^) for any ji € D(A“). 
Therefore, for t, s € [0, T], we have 

|j(e-*^-/)A“e-*^a:ll < 

< CpCa+pi^s-^^+^^xl (6.2.10) 
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We use (6.2.10) to obtain the inequality 

f ' II All ( |l/„(s. u)|| 

Jo 

■ [ l^(s -7’)| Il5n(7',«)ll dr] ds 
Jo 

f ' ||(e-<‘>-‘i)-< -/)e-("-)'>A“|| [ ||/„(s,u) 
Jo 

+ [ |a(s - r)| l|5n(r,u)l| dr] ds 
Jo 

< Ca,fi{t2 — tiY 


+ 


< 


( 6 . 2 . 11 ) 


where 


rpl-{<x->rS) 

From the inequalities (6.2.8), (6.2.9) and (6.2.11), it follows that {Snu){t) is contin- 
uous from [0, T\ into D{A°‘) with respect to the norm || ||a. Now, applying and 
subtracting from both the sides in (6.2.6), we get 


|l(S„K)(t) - ^)|U < ||(e-‘'‘-/)^|U+ ri|e-<‘->'‘A“|l[||/„(s,u)|| 

Jo 

+ [ |a(s-r)| ||5n('r,'u)|| dr] ds. (6.2.12) 

Jo 


It follows from (6.2.2-6.2.5) that 


< R. 


(6.2.13) 


Taking supremum over [0, T] we obtain that Sn maps Bii{Xa{T),^) into BR{Xa{T), ^). 
Now we show that Sn is a strict contraction on BR{Xa{T), ^)- For u,v E BR{Xa{T), ^), 
we have 


\\{SnU){t) - {SnV)m\a < f [ ll/n(s> 11 

Jo 

-f- [ la(s - r) 1 Ij^n (r, u) - ^„(r, v) H dr] ds 

Jo 

< (FRiTo) + lta|[iP(o,ro)ll<5filli«(0.ro)) 
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Tl-a ^ 

~ ^llvair) < -^{RFRiTo) + i21|al|iP(o,ro)llC?Rl|i,(o,ro)) 

< ^l|w-'y|U„(T)- (6.2.14) 

Here we have used (6.2.4) to obtain the last inequality. Taking supremum over [0, T], 
we get 

ll5nU — ‘S'„u||x^( 7 ’) < -||u — 'y|Ua(r)- 

Thus SnU is a strict contraction on BR{Xa{T),^). Therefore, there exists a unique 
Un € BR{XaiT)T$} such that SnU^ = Un- Clearly, satisfies (6.2.7). This com- 
pletes the proof of the proposition. 

Proposition 6.2 Let (HI) — (H4) hold. If 4> £ D{A°‘) for some 0 < o < 1, then 
Un{t) G D{A^) for all t G (0,r] where 0 < /? < 1. Furthermore, if <f> G. I){A) then 
Un{t) G D{A^) for all t € [0,T] where 0 < /0 < 1. 

Proof: Prom the Proposition 6.1 we have the existence of a unique G B^fX^ffT), 
satisfying (6.2.7). Part (a) of Theorem 6.13 in Pazy [21] on page 74 implies that 
g-tA . for t > 0 and 0 < ,5 < 1. Also, from Theorem 2.4 in Pazy [21] 

on page 4, we have € T>(A) if x G D{A). The results of the proposition follow 
form these facts and the fact that D{A) C D{A^) for 0 < < 1. 

Proposition 6.3 Let (HI) - (H4) hold. Then, for any <j> G T)(A“), 0 < a < 1 and 
any to G (0,T] there exists a constant Uta, independent of n, such that 

\\Unmp<Ut„ 0<^<1, to<t<T. 

Moreover, if <j> Q D{A), then there exists a constant Uq, independent on n, such that 

\\Unit)\\p < Uo, 0<^<1, 0<i<r. 

Proof: Applying A^ on both the sides in (6.2.7) and using part (c) of Theorem 6.13 
in Pazy [21] on page 74 after taking the norm, we have, for to ^ t ^T, 

rjil-p 


(6.2.15) 
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If (f>e D{A), then (f> € D{A^) for 0 < < 1 and we get 

ll«n(i)||^ < M\\A^<j>\\ + CfiL{R)^^!^<Uo. ( 6 . 2 . 16 ) 

This completes the proof of the proposition. 


6.3 Convergence 

In this section we prove the convergence of the solution u„ € Xa(T) of the approxi- 
mate integral equation 

Unit) = e-*^<t)+ [ 

Jo 

+ / a{s -T)gn{T,Un) dr] ds (6.3.1) 

Jo 

to a unique solution u of (6.2.1). 

Proposition 6.4 Let (HI) — (H4) hold. If <f) ^ L)iA°‘), 0 < a < I, then for any 
to € (0, T] 

lim sup l|u„(t) -1x^(4) lla = 0. 

{n>m, to<t<T} 

Proof: For n > m, we have 

l|/n('^) W?i) ~ /m(^! ^m)|| ^ ||/n(fj ^n) “ /n(^) ^m)!! 

+ WfnitjUm) ~ fm{tiUm)\\ 

< Fr{To)[ |lu„(t) - Um{t)\\a 
-h l|(F"-P’")u„^(f)|U]. (6.3.2) 

For 0 < a < < 1, we have 

< ( 8 ' 3 - 8 ) 

A77I 

Using (6.3.3) in (6.3.2), we get 

\\fnit,Un) - fmit,Um)\\ < Fk(Po)[ ll«^n(f) “ ^im(i)lU 

Am 


(6.3.4) 
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Similarly, we have 


Il5n(^> ~ 9miii %n 


Now, for 0 < io < io, we have 


\0-Oi ll^wWil/S]- 


(6.3.5) 


\Un{t) - Um{t)\\a < ^ [ |l/„(s, «„) - /,(5, n 

+ [ |o(5 - r)\ lj 5 „(r,'U„) - 5m ('r,Wm) II d'^] ds. (6.3.6) 
Jo 


The first integral can be estimated as 

rt'c\ 


Jo 

+ / |a(s - r)l Il5„(r, u^) - 5m (r, Mm) 11 dr] ds 
Jo 


< 2L{R)Ca{to-Q'%- 


For the second integral, we have 

f ||e-(‘-‘)'‘.4“|| 1 %) - /„(t, t<™)ll 

Jt'Q 

+ f |(2(s — t)| ||5^n('^)^n) ~ 

Jo 

Uif 

+ L{R)Ca _ g ja 11^” "" 

Using the inequalities (6.3.7) and (6.3.8) in (6.3.6), we obtain 

Iftin — '<Jm\\Xo,it) ^ 2L(R)CaitQ — ^o) ^0 

Ui> T^-° 

+ mc.^— 


a 


/** 1 

+ L(R)Ca J - ■'^mlU„(s) ds. 


(6.3.7) 


(6.3.8) 


(6.3.9) 
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Applying Gronwall’s inequality in (6.3.9) then letting m -)■ oo and using the fact 
that I'q is arbitrary, we get the required result. This completes the proof of the 
proposition. In view of Propositions 6.2 and 6.3, we have the following result. 

Corollary 6.1 //<^ e D{A), then 

lim sup |Mt)-Mi)IU = 0. 

{n>m,0<t<T} 

For the convergence of the solution u„(t) of the approximate equation (6.3.1) we 
have the following result. 

Theorem 6.2 Let (Hi) — (H4) hold and let (j) € D{A°‘). Then there exists a 
function u € Xa{T) such that Un u asn oo in X“(T) and u satisfies (6.2.1) on 
[0, T]. Furthermore, u can be extended to the maximal interval of existence (0, tmax)} 
0 < tmax ^ oo satisfying (6.2.1) on [0,tmaa:) o-^tl u is a unique solution to (6.2.1) on 
[0, tmax)’ 

Proof; We first assume that <j) € D{A). Corollary 6.1 implies that there exists 
u € Xa{T) such that u„ converges to u in A’q(T). Since Un € BR{Xa{T), $) for each 
n, u is also in Br{X^{T),$). Further, we have 

Wfnit, Un) - fit, «(i))ll = P^'^nit)) - fit, u(t))ll 

< FRiTo)[\\unit)-uit)\U 

+ |i(F”-/)u(t)|U]. (6.3.10) 

Taking supremum over [0,r], we get 

sup ll/n(ijUn) - /(t,u(i))j| 

0<t<T 

< Fh(To)[ i|u„ - ul|x„(r) 

+ ||P" - J|| ||ul|xo{T)] "^0 ss n -> oo. (6.3.11) 

Similarly, for g we have 

\\gnit, Un) - 9it, u(i))ll = IU(*> •f’”Un(t)) - git, u(t))|l 

< (?«(*)[ |lun(i)-u(i)lU 
+ jl(F"-/)u(i)|U]. 


(6.3.12) 
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From (6.3.12), we have 


o<t<r/ Il5n(s,«,i) - 5 ( 5 , 14 ( 8 )) II ds 

— ^i^)[ Ihn ~ '^i|Uo(r) 

+I1P” - /||M||x<,(r)] -> 0 as n -> 00 . (6.3.13) 

Using (6.3.10) and (6.3.12) in (6.3.1), we get 

u{t) = e-‘^0+ [ e-('-^)^[/(s,n(s)) 

Jo 

+ a{s -T)g{T,u{T)) dT]ds. (6.3.14) 

Now, let 0 c D{A°‘). Since for 0 < t < T, A°‘Un{t) converges to A°‘u{t) and 
Un(0) = 14(0) = <^, we have, for 0 < t < T, A“u„(t) converges to A^u{t) in H. 
Furthermore, since each Un is in BR{Xa{T),^) we have u e Br{Xc,{T),$) and for 
any 0 <to <T, 


lim sup ||wn(i) - w(i)lla = 0- 
"-^00 to<t<r 


(6.3.15) 


Also, 


and 


sup ||/„(t,'U„) - /(t,u(t))|| 

to<t<T 

< Fr{To)[ llun - 'u||x„(r) 

+ IjP" - J|| l|u!|xa(r)] ->-0 as n -)■ 00 (6.3.16) 

sup f |a(t-5)| ||5n(s,Un) “ 5 (s> u(s))|l ds 
to<t<T Jo 

< P(P)[ ||u„ - ut|xa(r) 

+ l|P"-lNk(r)]-^0 as n-4o<6.3.17) 


Now, for 0 < to < t, we may rewrite (6.3.1) as 

pto pt 

Ufi (t) = 6 




+ 


f a(s - r) 5 „(r, Un) dr] ds. 
Jo 


(6.3.18) 
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For the first integral, we have 


Thus, we have 


1/71(5, u„) + a{s-T)gn{T,Un) dT]ds\\ 

< ML{R)to. (6.3.19) 

\\Unit)-e-^U - re“^‘-*)^[/n(5,Un) 

JtQ 

+ / a(s - T)gn{r, u„) dr] ds|| 

< ML{R)to. (6.3.20) 

Letting n -4 cx) in (6.3.20), we get 

Jto 

+ a{s — T)g{T, u{t)) dr] ds|| 

^0 

< ML{R)tQ. (6.3.21) 

Since 0 < to <T is arbitrary, we obtain that u satisfies the integral equation (6.2.1). 
If u satisfies (6.2.1) on [0,T], then we show that u can be extended further. Since 
0 < To < 00 was arbitrary, we assume that 0 < T < To. We consider the equation 

dw{t) 


dt 


+ Aw{t) = F{t,w{t)) 


+ / a(t - s)G{s, w{s)) ds, 0 < t < To < oo, (6.3.22) 

Jo 

u;(0) = u{T), 

where F, G : [0, To - T] x T>(A“) H axe defined by 

F{t,x) = f{t + T,x) + h{t), 
rT 

h{t) = I o-it + T — s)g{$, u{s)) ds, 

Jo 

G{t,x) = g(t + T,x)j 



6.3 Convergence 


83 


for (f, ;J.) € [0,To — T] x D(/4“). We note that F and G satisfy 

(H3) and (H4), respectively, for To replaced by Tq — T. Hence there exists a 

^ ^ D{A )) for some 0 < Ti < Tq — T satisfying the integral equation 


w{t) = e ^\{T)+ [ 

Jo 

+ a{s- r)G(r, w{t)) dr] ds, 0<t<Ti. 
Jo 


(6.3.23) 


We define 


u{t) 


u{t), 0<t<T, 
w{t-T), T<t<Ti+T. 


Then u satisfies the integral equation 


5(t) = e-“^+ /‘e-(‘-)'‘[/(s,a(s)) 

Jo 

+ j a{s — T)g{r, u{t)) dr] ds, 0 <t <Ti+T. (6.3.24) 
Jo 

To see this, we need to verify (6.3.24) only on {T, Ti + T]. For t € [T, Ti 4- T], 


u{t) = w{t — T) 

= + / 
Jo 


t-T 


-(t-T-s 


'^^[F(s, w{s)) 


+ f a{s- t)G{t,w{t)) dr] ds 
Jo 


= e 


-it-T)A 


e-^^cl> + e-(^-)^[/(s,«(s)) 

+ f a,{s- T)g{T, u{t)) dr] ds 
Jo 

+ / a{s - t)G{t,w{t)) dT]ds 

Jo 





(6.3.25) 


Putting s + T = r] and r + T = ^, we have 

u{t) = / e"<*"*^^[/(s,«(s)) 

Jo 

+ [ a{s- T)gir, »{r)) dr] ds 

Jo 
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+ 


— g 


+ [ +h{ri - T) 

JT 

[ O’iv d^]dr] 

JT 

<j>+ re-(‘-*)^/(s,u(s))ds 
Jo 

+ / /* — r)g{T,u{T)) dr ds 

Jo Jo 

+ [ I — T)g{T,u{T)) dr ds 

Jt Jo 

+ f [ — T)g{T,u{T)) dr ds 

Jt Jt 

= e-‘U+ 

Jo 

4- / a(s — r)G(T, w(r)) dr] ds. 

Jo 


(6.3.26) 


Thus, u{t) satisfies (6.2.1) on [0,Ti +T]. Hence, we may extend u{t) to the maximal 
interval [0,tmax) satisfying (6.2.1) on [0,troax)j 0 < tmax < oo. 

Now we show the uniqueness. Let Ui and U 2 be two solutions of (6.2.1). Let T 
any number such that 0 < T < tmax- Let 

R = max{|l«i-(^||xa(T),||w2-^|U<,(r)}- 

Then, for 0 < r; < T, we have 

\\uiiv) - U2ir])\\^ < r i| [f{sMs)) - f{s,U2{s))\\ 

Jo 

+ [ |o(s - r)| l|^(r, Ui(r)) - g{T, U 2 (r))ll dr] ds 

Jo 

< L{R)Ca [ {r)- s)"“ll«i - “ 2 lUa(s) ds. (6.3.27) 
Jo 

Taking supremum on 0 < t? < t in (6.3.27), we obtain 

IK-«2|k(t) < L{R)C^ r{v-s)-%^i-U2\\x.is)ds. (6.3.28) 

Prom Gronwall’s inequality and the fact that 

llni(t)--W2(«)ll < - «2lk{T), 
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it follows that Uj — ^2 on [0, T]. Since 0 < T < tmax was arbitrary, we have Ui = «2 
on [O^ tffidx ). This completes the proof of the theorem. 

Corollary 6.3 If there exists a continuous function K : [0,oo) — > (0, oo) such that 

HOlU < K{t), t>0, 
then u satisfies (6.2.1) on [0, oo). 


6.4 Faedo-Galerkin Approximations 

For any 0 < T < tmax-> we have a unique u 6 X°‘{T) satisfying the integral equation 

u{t) = e~^^(t)+ [ e"^‘“®^^[/(s,u(s)) 

Jo 

+ f a{s - r)g{T,u{T)) dr] ds. (6.4.1) 

Jo 

Also, we have a unique solution € A“(T) of the approximate integral equation 

Unit) = e-^^<t>+ re-(‘-*)^[/(s,FX(s)) 

Jo 

+ f ais - T)g{T, P^Unir)) dr] ds. (6.4.2) 

Jo 

If we project equation (6.4.2) onto Hn, we get the Faedo-Galerkin approximation 
Unit) = P'^Unit) satisfying 

u[t) = / e"<*“*^^F"[/(s,«n(s)) 

Jo 

+ f ais- r)giT,Unir)) dr] ds. 

Jo 

The solutions u of (6.4.1) and u„ of (6.4.3), have the representations 

OO 

^i(t) = Y,a,it)ui, aiit) = iuit),Ui) i = 0,l,--- 


(6.4.3) 


(6.4.4) 


i =:0 

n 


Mt) = " («.(«), “i) i = 0,l.' 


(6.4.5) 
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Using (6.4,5) in (6.4.2), we obtain a system of first order integrodifferential equations 
da^(t) 

-^ + Aia"(t) = i?(i, «;(«).■■■, <W) 

+ / o(«-s)G?((, «;(«),■• ■,<(*))* (6.4.6) 

Jo 

a"(0) = (pi, (6.4.7) 


where 

i?{i,a5(t), ••.,<(«)) = lf(t,'£of{t)u,),u, 

\ i=0 

G“(i, <((),■ -,<(()) = 

and <pi = {<p,Ui) for i = 1,2, For the convergence of af to cq, we have the 
following convergence theorem. 



Theorem 6.4 Let (HI) - (H4) hold. Then we have the following, 
(a) If <p £ D{A°‘), then for any 0 < to <T, 


lim sup 


.1=0 


= 0 . 


(b) If <P£D{A), then 


lim sup 
n-^oo o<Kr 


5] A?“{ai(<) -<«)}’ 

Li=0 


= 0 . 


Proof; 


A‘-lu{t)-%,(t)] = A” 




,i=0 




i=0 


(6.4.8) 


(6.4.9) 


Thus „ 

II A“[u(t) - tin(t)]f > S 

i=0 

To conclude the results of the theorem, we may use the following result a proof of 
which is fiinni1a.r to the jKOofe of Proposition 6.4 mid CJorollaiy 6.1. 
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Proposition 6.5 Let (HI) - (H4) hold and let T be any number such that 0 < 
T < imax’ Then we have the following. 

(a) If <l> ^ then for any 0 < to <T, 


lim sup 

{n>m, to<t<T} 


1I^“M0-Wm(i)] II 


(b) If(j)eD{A), then 


= 0 . 


lim sup 

{n>m, 0<t<T] 


imut) - II 


= 0 . 


6.5 Applications 

We consider the following class of integrodifferential equations 
Ut-Au = fi{t,u,Vu) 

+ [ a{t- s)f2is,u{s),Vu{s)) ds, xGCl, t> 0 , ( 6 . 5 . 1 ) 

Jo 

u{t, x) = 0, X e do,, 
u{0,x) = tio(a^) 

Here f2 C R® is a bounded domain with sufficiently smooth boundary dQ and 
A, V are the 3-dimensional Laplacian and gradient, respectively. We assume that 
^ ^ and fi{t,u,p), {t,u,p) e R+ x R xR®, i = 1 , 2 , are locally Lipschit2 

continuous functions of ail its arguments. We further assume that there are fiinctions 
Pi : R+ X R4. R+ and a real constant 7> 1 < 7 < 3 sudi that 

\fi{t,u,p)\ < pi(tN)(i + bn. 

\fiit,u,p) - fi{t,u,q)\ < Pi(t,iw|)(l + br“^ + b|'^“^)|p-g|, 

\fiit,u,p) ~ fi{t,v,p)\ < Pi(i,N + M)(i + br)b-^l 

for all t for i = 1 and a.e. i for i = 2 where pi(.,r) is nondecreasing and p2(., r) € 
00) for each r > 0 where ( 1 /p) + (I/9) = 1 - 


We reformulate (6.5.1) as an abstract integrodifferential equation ( 6 . 1 . 1 ) in the 
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Hilbert space H = where A = —A + cl, for some c > 0 with D{A) = 

H2(f2) n Hj(fi), are Sobolev spaces (cf. §7.1 in [21] for the defini- 

tions) and I is the identity operator. 

We define the maps f,g : [0,oo) x D(A°‘) by 

f(t,u){x) = fi{t,u{x),Vu{x)) + cu(x), 
g(t,u)(x) = f 2 (t,u(x),Vu(x)), 

where max{3/4, (5j - 3)/47} <Q;<land(^eHis given by 4>(x) = uo(x}. We 
note (Hl)-(H4) are satisfied (for (H3) and (H4) cf. Theorem 4.4 pp. 244-245 
in Pazy [21]). Thus, we may apply the results of the earlier sections to guarantee 
the existence of Faedo-Galerkin approximations and their convergence to the unique 
solution of (6.5.1). 



Future Scope 


In this chapter we mention some of the problems which we plan to study in the near 
future. 

1. Evolution Systems and Applications 

All the problems we have considered in the present work, the operators ap- 
pearing in the abstract Cauchy problems does not depend on time. In many 
practical applications, the coefficients of the differential operators are in gen- 
eral depending on time. Therefore in the abstract formulations of such prob- 
lems, the associated operators will also be time dependent. We plan to study 
such problems with the help of some modifications in the ideas and techniques 
developed in the present work. 

2. Quasi-linear Evolution Equations 

More generally, as mentioned above, the coefficients of the partial differential 
operators will depend on the unknown solution of the problem. Therefore 
the associated operator in the abstract Cauchy problem will involve the un- 
known function. Many of such problems are in the domain of quasi-linear 
evolution equations. We plan to modify the techniques of the method of semi- 
discretization in time to study these quasi-linear evolution equations. 

3. Approximation Results for Hyperbolic Problems 

In chapter Six, we have considered the Fadeo-Galerkin approximation of solu- 
tions to an integrodifferential equation. This integrodifferential equation is an 
abstract formulation of a parabolic integrodifferential equation and therefore 
the assumption that the associated operator is the infinitesimal generator of 
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analytic semigroup. As mentioned earlier, this allows us to consider the frac- 
tional powers of the operator. For a hyperbolic problem, if we consider its 
abstract formulation as a first order equation in a Banach space, the associ- 
ated operator is the infinitesimal generator of a strongly continuous semigroup 
only. This restricts the analysis of such problems considerably. We plan to 
consider the abstract formulations of such problems as second order equation 
and the associated operator will still be generating analytic semigroup. We 
hope to modify our techniques suitably in order to study such problems and 
prove some approximation results. 
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